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HAUPTAUFSÄTZE 


Iterationsverfahren für komplexe Nullstellen algebraischer 
Gleichungen 


Von L. Collatz in Hannover 


Es wird ein Iterationsverfahren beschrieben, welches bei einer beliebigen Wurzel einer algebraischen 
Gleichung mit reellen oder komplexen Koeffizienten konvergiert, sofern man über einige, noch frei wählbare 
Größen passend verfügt und welches häufig rascher als das Newtonsche Verfahren Nüherungswerte für 
komplexe Nullstellen liefert, oft sogar für mehrere komplexe Nullstellen zugleich. 


A method of iteration is described that converges for any root of an algebraic equation with real or complex 
.coefficients provided that certain quantities, the values of which can be chosen freely, are given appropriate 
values. Frequently, Ihe method gives approximate values of complex zeroes more rapidly than the method 
of Newion, somelimes even several complex zeroes at the same time. 


Une möthode & it£ralion est donnee,,qui converge pour un zero quelcongue d’une Equation algebrique & 
coefficients reels ou complexes, si certaines quantites encore disponibles sont choisies convenablement. Souvent 
la methode nous fournit de zEros approwimatifs plus vite que la methode de Newton; frequemment elle möme 
donne plusieurs zeros complewes a la fois. 

OnmeziBaetcat METON HOCHENOBATENBHBIX TIPHONMKeHNa, CXonmmmich B cAyyae 10Doro 
KOPHA aureöpanyeckoro YPABHeHHA C BEINECTBEHHEIMH HIIH KOMINIEKC HEIMU KOSPHLMEHTAMH, 
eCJIH HEKOTOPEIE HOCTOAHHBIE, BEIIHYHHA KOTOPEIX MOKeT ÖBITb IPOHSBONBHO NPuHATa, OyayrT 
HOANYOANSIIEe TONOÖPAHBI. HTOT METOA AAeT yacTo ÖkIcTpee Merona HBWTOHA IPHÖNNSKeHHBIE 
3HAYeHHA WI KOMIITEKCHBIX KOPHeH, HEPENKO AAHce MIA HECKONBEHX KOMILTEECHEIX KOPHei 
OAHO BpeMeHHo. 


An Näherungsverfahren zur Bestimmung der Nullstellen algebraischer Gleichungen be- 
steht kein Mangel!). Insbesondere läßt sich, wenn ein Näherungswert für eine Wurzel bekannt 
ist, mit dem Newtonschen Verfahren unter Verwendung des Hornerschen Schemas auch 

bei komplexen Wurzeln die Genauigkeit leicht bis zu einem alle technischen Bedürfnisse befrie- 

 digenden Grade steigern. Mitunter aber macht es Mühe, besonders bei Vorhandensein mehrerer 
Paare komplexer Wurzeln oder bei komplexen Koeffizienten in der vorgegebenen Gleichung, ge- 
eignete Näherungswerte zu erhalten, und dabei kann, wie die folgenden Beispiele zeigen, auch das 
Newtonsche Verfahren sich sehr ungünstig verhalten; daher sei ein Verfahren beschrieben, 

welches in vielen Fällen rascher als das Newtonsche Verfahren zum Ziel führt und dabei 
häufig zugleich für mehrere Nullstellen Näherungen liefert. 


1. Das allgemeine Verfahren 
Vorgelegt sei die algebraische Gleichung 


I) Sun" tn "+... 49240, = 0 A, SER Berl oh s (1) 
e= 


0 mit zahlenmäßig gegebenen reellen oder komplexen Koeffizienten a,. Es sei +0, +0; die 
- Wurzeln £, sind gesucht. Das Iterationsverfahren besteht darin, daß man die Gleichung mit 
{. einer Potenz 2” multipliziert, dann nach 2? oder allgemeiner nach (z* — q)? auflöst, also auf 
_ die Gestalt bringt ; 
ee N... (2), 


_ wobei R(z) dann eine bekannte rationale Funktion von z ist, und sodann von einer gewählten 
Zahl 2, ausgehend eine Zahlenfolge z, bestimmt, für welche 


N 2; (@F,,— N’ =R@%) N EL El) 


gilt. Dabei wird im allgemeinen z,,2,,... mehrdeutig, und zwar z, wird (kp)‘.Werte annehmen 
_ können. Trägt man sich die z, in der komplexen Zahlenebene auf, so erhält man eine Punkt- 
menge, die im allgemeinen außer den ö, noch weitere Häufungspunkte haben wird, bei der man 


Mn — 


z.B. Fr. A. Willers: Methoden der praktischen Analysis. Berlin und Leipzig 1928. Dort sind 
lere S. 194ff. auch ausführliche Methoden zur Berechnung komplexer Nullstellen behandelt. 
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aber normalerweise schon beim Aufzeichnen der Punkte z, und 2, sehen wird, welche der (kp)? 
Werte z, als Näherungen für eine Wurzel in Betracht kommen, vgl. die Beispiele. Dieses 
Iterationsverfahren enthält einige der bisher benutzten als Spezialfälle. 


2. Spezialfall des Reellen 
Im Reellen lassen sich die Verhältnisse leicht überblicken daher sei es gestattet, zunäehst 
den Fall, daß alle a, und eine der anzunähernden Wurzeln £, reell sind, herauszugreifen. 
Für diese Wurzel schreiben wir kurz &. Wegen a,+0 ist &+0. Die Gleichung (1) werde mit 
z’ multipliziert und sodann nach einem Gliede mit dem Koeffizienten a„, der etwa +0 sei, 
aufgelöst: 
N 
ENT Er un.) Pe 
v=0 
v=m 
Der Strich an dem Summenzeichen soll dabei auf das Auslassen des Gliedes mtv = m 
aus der Summe hinweisen. Die Iterationsvorschrift (3) soll dann lauten 


„ir =) a De ae us BER 2 05h; 


mit 
1 e ER 
(2) a NUT Eh Te Be re On 
Pe) | a ( 
v==m 


Über die eventuelle Mehrdeutigkeit der Vorschrift muß eine besondere Festsetzung ge- 
troffen werden, worauf wir noch zu sprechen kommen. 

Bekanntlich konvergiert ein solches Verfahren in der Nähe der Wurzel, wenn !p'(&)|< 1 
ist, und zwar konvergiert es um so besser, je kleiner |p’(&)| ist. Die Ausrechnung ergibt wegen 


ol)=L 


—1 mr % v+r—1 en! n { 
! \ 5 y‚r . N -—m yr a v 
P(I)= Bar =v +rn)aÖ See mo tr) 
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Wegen Fl) Da, r=0 
0 
folgt (am Summenzeichen bleibt jetzt der Strich fort) 
; En n BAR, | n 
yi)—= mi +r)a,L& a PT N gr 


erde (m-+-r) v=0 — Amim + rn)... »v=0 
oder schließlich 


ed)=1 a -(M) 

mit 
Sa 
Cd) a mai . (8). 


‚ Nun setzen wir die reelle Zahl C„, (£)+0 voraus; ist nämlich (m (&) = 0, so hat f(z) = 
bei 2=[ eine Doppelwurzel, und es läßt sich nach den Methoden der Algebra eine Gleichung 
niedrigeren Grades herstellen, die z={£ als einfache Nullstelle hat, für die also die entsprechende 
Größe C„(2) von Null verschieden ist. Für C,(£) +0 läßt sich die Größe r so wählen, daß in (7) 
\P(&)l< 1 wird und das Iterationsverfahren konvergiert. Es gibt dann ein günstigstes r, so daß 
o'(&)=0 wird, nämlich 

Fate) 


a BR 


. Hat man Näherungswerte für £ und f’(£), so kann man auch r geeignet wählen, normaler- 
weise wird man Topt. durch eine ganze Zahl annähern. ‚ 

Durch die passende Wahl von m und r (und durch etwaige Abspaltung mehrfacher Null- 
stellen) läßt sich also bei einer reellen Nullstelle einer algebraischen Gleichung mit reellen Ko- 
effizienten stets erreichen, daß das Iterationsverfahren (5), (6) konvergiert. 


r=rp =— m— („(t)=—m+ 


3. Konvergenzbeweis bei komplexen Wurzeln a 
Die Überlegungen von Nr. 2 lassen sich ins Komplexe übertragen. Es seien also die a,und 


_ 


- Z beliebig komplex, jedoch sei Z eine einfache Nullstelle, also f’(£)+0, was sich nach den obigen 
‚Bemerkungen stets erreichen läßt. r soll aber reell bleiben. Die Formeln (5), (6), (7), (8) gelten 
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auch jetzt, und wieder ist eä)l< 1 hinreichend für die Konvergenz des Iterationsverfahrens, 
wenn der Ausgangswert z, genügend nahe bei £ liegt?). 

Die rechte Seite von (6) ist eine mehrdeutige Funktion. In einem bestimmten der (m + r) 
Blätter der zugehörigen, geeignet aufgeschnittenen Riemannschen Fläche wird für 2={& 
wieder der Wert Ö angenommen. Die Stelle 2 = Zist wegen {+0 kein Verzweigungspunkt. Wir 
denken uns die (m =E r)-te Wurzel so gezogen, daß man in diesem bestimmten Blatt der Rie- 
mannschen Fläche bleibt. 

Nun denken wir uns für alleWerte von m, für die «m +0 ist, die zugehörigen Zahlen C (£) 


nach (8) berechnet und die Zahlen 


CO (£)#0). Nun ist 


in der komplexen Zahlenebene dargestellt (es ist 


1 
Cm (2) 
also 


wobei der Doppelstrich am Summenzeichen bedeuten soll, daß nur über die Werte » zu sum- 
mieren ist, für die a, von Null verschieden ist. Es können, da die linke Seite der Gleichung reell 
ist, insbesondere nicht alle C',(£) rein imaginär sein. Betrachtet man nun die Werte von 9 (£) 
n (7) für verschiedene reelle Werte von r, so liegen die zugehörigen Punkte 9’ (£) in der komplexen 
Zahlenebene auf einer Geraden, die durch den Punkt 1 geht, die aber für ein passendes m, für 
welches C'„ (2) nicht rein imaginär ausfällt, in den Einheitskreis eindringt. Für gewisse r liegt 
daher sicher ° (£) innerhalb des Einheitskreises. Somit gilt: 


Bei einer beliebigen komplexen (als einfach voraussetzbaren) Wurzel einer algebraischen Glei- 
chung läßt sich a durch passende Wahl von m und einer reellen Zahl r und durch passende Fest- 
legung über die in (6) auftretende Mehrdeutigkeit erreichen, daß das Iterationsverfahren (5), (6) kon- 
vergiert, sofern der ee z, der Wurzel £ genügend benachbart ist. 


4. Beispiele 
Zunächst seien einige allgemeine Bemerkungen über den Vergleich mit dm Newton- 
schen Verfahren gemacht. 
Wenn man das Newtonsche Verfahren anwendet 


Zm+1 — m a) 


und von einem Wert z, ausgeht, der von allen Wurzeln £, weit entfernt ist, so ist überschlags- 
weise 


fi) =mII (2, — Ce) 7 An? , 


0= 
I) zn Ar 


J (zo) w2 el 
W 


» 
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1 
. d.h. der neue Wert z, liegt grob gesprochen um 55 der Entfernung näher an den Wurzeln. Bei 


größeren Werten von n bedeutet das eine sehr langsame Konvergenz. Durch Wurzelziehen 
kommt man im allgemeinen schneller in die Nähe der Nullstellen. In den folgenden Beispielen, 
in denen das hier beschriebene Iterationsverfahren mit dem New tonschen verglichen wird, be- 
stätigt sich das durchaus. 
I. Gleichung dritten Grades: 
—i2+2=0. 
Hier lösen wir etwa nach dem al mit 2? auf und erhalten die Iterationsvorschrift 
Ara = Yi 2, —2. 
Das Iterationsverfahren läßt sich hier sehr bequem zeichnerisch durchführen. Es werde z. B. 


a gewählt, und dann die drei Werte z, ermittelt. Bei der Aufstellung von 2, zeigt sich, daß 
- jeweils einer der drei neuen Werte z, in der Nähe des zugehörigen Ausgangswertes 2, liegt, ent- 


2) Vgl. Fr. A. Willers: a2. 0. 5. 165. 
ie 


Zi 
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hand bei z,, so daß die Festlegungen der dri | 
hat also hier, vgl. Abb. 1, gleich alle drei Wurzeln erhalten. Die Zeichnung 


a 0,85 +1,21 3 


Er & =! — 1,26—0,27 3 


un" -0,41-0,95 3 
Zaof >% 


j@)=?—-i2+2=0 & 

8 z 
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Eine Verbesserung nach dem Hornerschen Schema liefert dann die genaueren Werte 


0,8561 4'1,22474 mn Rn N U ar Se 
& = — 1,2608 — 0,2684 i DR NER 
0,4046 — 0,9562 8... una Le 
In Abb.1 sind ferner die nach dem Re a ae 
Newtonschen Verfahren erhaltenen N vo Bo) 2 > 


Werte z, bis z, eingetragen, die man von 
demselben Ausgangswert z, = 4 aus er- 
‚hält; die Abbildung zeigt, daß bis zum 


Br: je=z#—2ie +2. 3:i=0 a Le: 


MT Iterationsverfahren 


2,41 = 3li 3 Yo, +84—1 
— "—  — Newtonsches Verfahren 


2,) 
I een; 


Ef "© Nullstellen ‘ 
«@=sH+iy) i 
Abb. 2 
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Wert 2; noch von keiner Annäherung an irgendeine der Wurzeln gesprochen werden kann. 
Man wandert bei Anwendung des Newtonschen Verfahrens gemächlich zwischen den 
Nullstellen hindurch. Auch wenn man von einem anderen Wert, etwa z, = 1 ausgeht, zeigt sich 
zwar eine Konvergenz gegen eine Wurzel, aber der hier erhaltene Wert z, ist noch viel schlechter 
als der Wert z, nach dem hier beschriebenen Iterationsverfahren. Das hängt natürlich davon ab, 
wie gut der Anfangswert z2,ist. Dieses Beispiel soll nur zeigen, daß das Newtonsche Verfah- 
ren ungünstig sein kann, wenn man keine geeignete Näherung kennt. 


II. Man kann durch eine Transformation 2 = 2° — en in (1) das Glied mit der (n—1)- 
n*a 


- Rn 
ten Potenz von 2’ beseitigen und dann die Glieder mit 2’* und 2’*-2 zu einem vollständigen 
Bat ergänzen, wodurch man häufig bei der Iteration (3) bessere Ergebnisse erhält. 
ei | 
#—212 +2 —3i=0 

schreiben wir 

(?— 1? +2—3i+1=0 
und iterieren nach der Vorschrift 


= til mt —l. 
Die Durchführung für den Anfangspunkt 2, =1 ist in Abb. 2 dargestellt. Bereits bei 2, 


hat man vier Wurzeln gut angenähert, und die Überlegenheit gegenüber dm Newtonschen 
Verfahren zeigt sich hier ebenso wie in Abb. 3 bei dem Beispiel einer Gleichung 5. Grades. 


Eingegangen am 25. Mai 1949. 


Näherungstheorie des unsymmetrischen Schalldurchgangs 


in einer Lavaldüse 
Von Hermann Behrbohm, Waldkirch/Brsg. 


Das Verhalten einer ebenen Düsenströmung in der Umgebung des kritischen Zustandes wird approwi- 
mativ durch die Differentialgleichung der Wärmeleitung beschrieben und deren Lösung bei den entsprechen- 
den Rand- und Anfangsbedingungen angegeben. Durch geeignete Spezialisierung der Einlaufsbedingungen 
wird ein einfaches Beispiel erhalten, das alles physikalisch Wesentliche bereits zeigt. 


_ The author gives an approximate description of a plane flow in a nozzle near the critic state by the diffe- 
renlial equation of heat iransfer and shows the solution for adequate initial and boundary conditions. By 
an appropriate specialization of the entrance conditions, he obtains a simple example showing all that is 
physically interesting. d z 

L’auteur decrit d’une facon approchee le comportement d’un &coulement plan dans une tuyere en prowi- 
mite de V’e&tat crilique par l’öquation differentielle de la chaleur et indique la solution pour les conditions 
initiales et aux limites correspondantes. On obtient'un ewemple simple par une sp£cialisatiyn convenable 
des cynditions d’entröe qui met en Evidence toul ce qui est d’une importance physique. 

Ilpu nomoımu auhepeHumasısHoro YyPaBHeHHA TEILMIONPOBOAHOCTH IPHÖNSKEHHO OTHMCHIBACTCH 
INIOCKO® TEYEHHE B BEIITYCKHOM OTBEPCTHH IPH IIPHÖJNHSKeHHU K KPHTHYeCKOMY COCTAHHIO MH 
AaeTca pemeHue MIA COOTBETCTBEHHBIX TPAHHYHBIX U HAYAJIbHBIX YCAOBHH. Cmenmarusupys 
HOAXONHINUM O0Pa30M BXONHBIE YCIOBHA, ABTOP HONYYAeT IPOCToä HPHMePp, TORAsEIBamımma 
Bce& Pusnyecku CYIH&ECTBEHHO®, 


1. Das Verhalten einer ebenen reibungs- und wirbelfreien Düsenströmung in der Um- 
gebung des kritischen Zustandes ist bereits verschiedentlich Gegenstand der Untersuchung 
gewesen [1], [2], [3], [4]. Die befolgte Methode war dabei stets die des Potenzreihenansatzes 
in der Strömungsebene mit anschließender Ermittlung endlich vieler Anfangskoeffizienten der 
Reihen. Untersuchungen über die Konvergenz der verwendeten Reihen wurden nicht durch- 
geführt. Im folgenden wird ein auf wesentlich anderen Prinzipien beruhendes Näherungsver- 
fahren zur Beschreibung des Strömungsfeldes in der Umgebung der engsten Stelle einer Laval- 
düse bei maximalem Durchfluß angegeben). 

Man lege dazu die ebene, reibungs- und wirbelfreie Strömung eines adiabatisch kompres- 
_  siblen Mediums ‘zugrunde. &, y seien die rechtwinkligen Cartesischen Koordinaten in der 
_  Strömungsebene, v, v die Komponenten der Geschwindigkeit w = Yu® + vo in Richtung der 
 &- bzw. y-Achse.. Es seien ferner p und o der Druck bzw. die Dichte des strömenden Mediums. 
Die örtliche Schallgeschwindigkeit werde mit ce bezeichnet, so daß c gemäß der Adiabate durch 


#5 definiert ist (2 = 0)/% = Verhältnis der spezifischen Wärmen bei konstantem 


i 1) Herrn K. Oswatitsch danke ich auch an dieser Stelle für anregende Diskussionen bei der Nieder- 


dieser Arbeit. ö 
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Wird hierin (mit M= 
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Druck und konstantem Volumen). Der durch = « Pesehene Zustand heißt der ie he i 
Zustand, alle im kritischen Zustand betrachteten Größen heißen ‚‚kritische‘“ Größen. Sie werden A 
durch einen Stern bezeichnet. 

Die Forderung der Wirbelfreiheit und die ern der Stetigkeit den NE SFURTRER, er- 
lauben die Einführung von Geschwindigkeitspotential (x, y) und Stromiuplsiiitir (2, Y), die x 


bis auf jeweils eine additive Konstante durchess . Ki 
U [ ® 203 5 ‚ 
x == er ; : Py er A ee NE: } . WR (1) Ä 
bzw. k: 
e® a .. 
Ye —— pre FE. | | Yy = o* c* a > res (2) 
definiert sind. Damit wird | 
3(9, Y) [2 Aa ar 
= — een uals.... 6 N RO 
ER = a — Pte FANG; 8). 


also D ++ 0 für ow2 +0. Von etwaigen Staupunkten und vom Vakuum abgesehen, ist daher das 
durch = o(z, y), y = y(«, y) vermittelte Bild der z, y-Ebene auf die 9, y-Ebene eineindeutig, 
so daß die Strömung auch ebenso gut in der 9, y-Ebene wie in der x, y-Ebene studiert werden kann. 

Neben dem Geschwindigkeitsbetrag w werde noch der Strömungswinkel © = arctg v/u ein- 
geführt. Der Zusammenhang zwischen «, y-Ebene einerseits und @, y- us Reel wird 
dann durch die Differentialbeziehungen 


f} 


d2=°. 00s0dp X sin@dy v| 
De 
Ä . (4) 
RR ER R\E | { 
Ay r aReT v 


gegeben. ’ ea 
Mit o und yals unabhängigen Veränderlichen lauten die Gleichungen der Drehungsfreiheit 7 
und der Kontinuität . 


2.2 (Or 
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C 


— Machscher Zahl) 


definierte Funktion w(; 2) eingeführt, so erhält man aus. o) 


iterentialkieichungen 
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“ für die beiden Funktionen W=W(p, y und 0 = cr 


folgenden Unterskenisfen dienen Boll; Der Faktor von — ist 
‘9 von el allein: | 
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Da das Hauptaugenmerk dieser Note sich auf die Umgebung des kritischen Zustandes richten 
soll, werde die gemäß (6) noch freie additive Konstante so gewählt, daß W(1) = 0 ist, d.h. 


es werde 
wic* ' w 
0 e 


W= (9) 


gesetzt. W < 0 bedeutet dann Unterschall-, W > 0 dagegen Überschallgeschwindigkeit. Wegen 


1 1 
0 FE s—1 ll. 
DIA NR #+1 \c* 


Zahlentafel 1 


w 


W (2) und K (2) nach Formel (1,10) 


w|c* W 14 | w/c* W K wjc* W K 

) es —0,40188 | 0,855 | —0,1679 | —-0,23570 1,11 0,0988 0,37044 
0,1 —2,9617 | —0,40186 | 0,56 | 0,1615 | —-0,23043 1,115 0,1029 0,39367 
0,15 | 2,6916 | —0,40179 | 0,865 | -:0,1551 | —-0,22500 1,12 0,1066 0,41761 
0,2 — 2,2435 | —0,40160 | 0,87 | 0,1482 | --0,21941 1,125 0,1105 0,44228 
0,25 | 1,8988 | —0,40119 | 0,875 | 0,1418 | —-0,21366 1,13 0,1143 0,46770 
0,3 —1,6214 | —0,40042 | 0,88 | 0,1356 | --0,20774 1,135 0,1182 0,49388 
0,32 | —1,5220 | —0,39998 | 0,885 | —-0,1292 | —-0,20165 1,14 0,1219 0,52088 
0,34. | —1,4310 | —0,39943 | 0,89 | 0,1229 | —-0,19538 1,145 0,1256 0,54869 
0,36 | —1,3454 | —0,39876 | 0,895 | --0,1159 | —-0,18893 1,15 0,1294 0,57736 
0,38 | 1,2649 | —0,39796 | 0,9 —0,1107 | —-0,18229 1,155 0,1331 0,60691 
0,40 | —1,1898 | —0,39701 | :0,905 -| —0,1044 | —-0,17545 1,16 0,1366 0,63737 
042 |—1,1173 | —0,39589 | 0,91 | 0,0990 | —-0,16841 1,165 0,1402 0,66877 
0,44 | —1,0495 | —0,39456 | 0,915 | —0,0925 | —-0,16117 1,17 0,1437 0,70114 
0,46 | —0,9851 | —0,39300 | 0,92 .| —0,0875 | —-0,15372 1,175 0,1473 - 0,73451 
0,48 | 0,9240 | —0,39119 | 0,925 | 0,0809 | —-0,14605 1,18 0,1507 0,76892 
0,5 —0,8659 | —0,38909 | 0,93 | 0,0751 | —-0,13815 1,185 0,1540 0,80440 
0,52 | 0,8105 | —0,38667 | 0,935 | —-0,0694 | —-0,13003 1,19 0,1576 0,84099 
0,54 | 0,7578 | —0,38388 | 0,94 | —-0,0637 | —-0,12167 1,195 0,1610 0,87871 
0,56 | —0,7048 | «0,3068 | 0,945 | 0,0582 | —-0,11306 1,2 0,1642 0,91762 
0,58 | 0,6591 | —0,37701 | 0,95 | --0,0525 |: —-0,10420 1,21 0,1707 0,99916 
0,60 | —0,6134 | —0,37283 | 0,955 | 0,0471 | —-0,095080 1,22 0,1772 1,0859 
0,61 | —0,5917 | —0,37052 | 0,96 | 0,0416 | —-0,085695 1,23 0,1834 1,1783 
0,62 | 0,5689 | —0,36807 | 0,965 | 0,0355 | —-0,076035 1,24 0,1896 1,2766 
0,653 | —0,5478 | —0,36544 | 0,97 | 0,0310 | -—-0,066092 1,25 0,1957 1,3813 
0,64 | —0,5267 | —0,36265 | 0,975 | —0,0257 | -—-0,055858 1,26 0,2017 1,4929 
0,65 | 0,5065 | —0,35968 | 0,98 -| —0,0205 | -—-0,045323 1,27 0,2073 1,6119 
0,66 | 0,4863 | —-0,35652 | 0,985 | —-0,0152 | -—-0,034479 1,28 0,2129 1,7387 
0,67 | 0,4669 | —0,35315 | 0,99 | 0,0108 | --0,023318 1,29 0,2184 1,8740 
0,68 | 0,4475 | —0,34957 | 0,995 | —-0,0049 | —-0,011828 1,30 0,2238 2,0183 
0,69 | 0,4292 | —0,34576 | 1,0 0 f) 1,32 0,2339 2,3370 
0,70 | 0,4108 | —0,34172 | 1,005 0,0050 0,012176 1,34 0,2440 2,7007 
0,71° :| 0,3930 | —-0,33741 | 1,01 0,0099 0,024710 1,36 0,2584 3,1166 
0,72 | 0,3751 | —0,33284 | 1,015 0,0148 0,037614 1,38 0,2625 3,5930 
0,73 | —0,3577 | —0,32800 | 1302 0,0196 0,050898 1,4° 0,2710 4,1399 
0,74 | —0,3403 | —0,32285 | 1,025 0,0243 0,064575 1,42 0,2792 4,7688 
0,755 | 0,3239 | —0,31738 | 1,03 0,0290 0,078656 1,44 0,2870 5,4938 
0,76 | 0,3074 | —0,31158 | 1,035 0,0340 0,093153 1,46 0,2944 6,3315 
0,77 | —0,2917 | —0,30543 | 1,04 0,0384 0,10808 1,48 0,3015 7,3016 
0,78 | 0,2760 | —0,29891 | 1,045 0,0446 0,12345 1,5 0,3082 8,4280 
0,79 | —0,2607 | —0,29199 | 1,05 0,0475 0,13928 1,6 0,3365 17,651 
0,850 | 0,2455 | —0,28465 | 1,055 0,0520 0,15557 Wr 0,3577 39,166 
0,805 | 0,2381 | —-0,28082 | 1,06 0,0566 0,17236 1,8 0,3728 95,016 
0,81 | 0,2309 | —0,27688 | 1,065 0,0610 0,18964 1,9 0,3832 262,58 
0,815. | —0,2236 | —0,27281 | 1,07 0,0652 0,20744 2,0 0,3898 878,91 
0,82 | --0,2164 | —0,26863 | 1,075 0,0696 0,22576 2,1 0,3936 | 3957,1 
0,825 | 0,2088 | —0,26432 | 1,08 0,0738 0,24464 2,2 0,3955 | 29552 
0,83 | 0,2013 | —0,25989 | 1,085 0,0781 0,26409 2,3 0,3961 | 62793.101 
0,835 | —0,1950 | —-0,25533 | 1,09: 0,0822 0,28411 2,4 0,3963 | 46702.104 
0,834 | 0,1883 | —-0,25063 | 1,095 0,0863 0,30474 v6 0,3963 4 
0,845 | —0,1814 | —0,24580 | 1,1 0,0905 0,32599 e 

0,85 1,105 0,0946 0,34788 


—0,1742 | —0,24082 


MERN 


ei 104 Behrbohm, Näherungstheorie des u J 1 
ie. | u A. 7 = 


# wird mit x = 1,4 für Luft als das hier interessierende Medium 
= w\2 1/w\alı/ If: 
W 0,7887 ) 2) Be. 
A | I -26 u) o* ı 6 \c 
K \ 8 
= 
K C 
er | 1/ w\2\6 
5 wit 
z14 Ll 


Das Sr (7) en für M< % d. ne W<o0 
von elliptischem, für M >1, d.h. für W>0 von 
hyperbolischem Typus. Das Interesse liegt auf der | 
beiderseitigen Umgebung von W=0. Wegen \ u 

} nn 
KO) 0 RN (5 =x%+1 
beginnt die Taylorsche de von K(W) in der 
Umgebung von W = 0 mit N 
E(M = (+1) W+:-- 
Ersetzt man also den Koeffizienten X(W) durch seine 


Tangente im Nullpunkt, so geht das System (7)inder 
Umgebung des Schalldurchgangs in das nichtlineare 


SEE 


Näherungssystem 
BEIN oW 
op  ay 
as br IR 

9 We Re 

uw | 

9Y oo 

w(4 as =) gemäß (1,10) für 0,05 4 S14 N 

ge, ns über. Es hat wie das Ausgangssystem (7) die Eigen- 
ER schaft, für W<0, d.h. im Unterschall, elliptischh 
B; für w >0 dagegen, d.h. im Überschall, hyperbolisch zu sein, und dient der zu entwickelnden a 


Näherungstheorie als Grundlage. 


2. Die erste Gleichung dieses Systems erlaubt die. Einführung eines allgemeinen Potentiale ' rs 
© == w(p,y) derart, daß W und © durch er 


= _..0@ _. 0@ 
a in PR TE Ir 1. nn ER RR Sekt = (1) 


gegeben werden. Gemäß der zweiten Gleichung genügt dann o(p, y) der nichtlinearen partiellen a 4 
Differentialgleichung 2. Ordnung. Ä * & 


sionalen  iine & De der sogenannten Siromtadenthear 1 
keit W in erster Näherung längs jeder Potentiallinienstrecke zwische i 
K y= y, als konstant, d. h. von 9 a Be a er 


als bekannte Funktion von @ sn angesehen: 


oW aW_ 
We 
Win) 9 
aribr man dann die zweite Gleichung (1, 1) nach % vo 


hiermit 
t Y 


00 | a" IE 
[ar = 9,W- 9 = +1 W 


MEN 
bi 
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oder also für den Gradienten @(p) der Stromfadengeschwindigkeit W (9) 


1 9,(9) — 9, () 
@ res a Be en N N EN . 
(m) Yı 9% (#+1)W(p) (3) 


Ersetzt man alsdann o,, - in (2) durch diesen Näherungswert (3), so erhält man die 
Näherungsgleichung ? 

Beh), GO en. en ine (4), 
so daß damit die nichtlineare Differentialgleichung (2) von gemischtem Typ eine lineare Differen- 
tialgleichung von durchweg parabolischem Typ geworden ist. 

Da sich das Interesse (der Natur der gewählten Näherung gemäß) im wesentlichen auf das 
Strömungsverhalten in einer Lavaldüse in Schallnähe konzentrieren soll, so soll y=y, als 
Düsenachse angesehen und auf ihr © = ©,(p) = 0 angenommen werden. Die Stromlinie y= y, 
dagegen entspreche der Düsenkontur, die dann durch die Wahl von W(p) wegen der aus (3) 
fließenden Beziehung 


9,(P) = (x +1) (y, — Yo) W(p) ae NE (5) 


bereits festgelegt ist, oder wenn man — je nach Standpunkt — so will, die nach ihrer Vorgabe 
in der Gestalt © = ©,(p) ihrerseits gemäß (5) den Gradienten @(g) in der Gleichung (4) festlegt. 
Die engste Stelle der Düse werde nach = p* >0 gelegt, so daß bei maximalem Durchfluß 


W(9*) =0 gilt. Die Stromfunktion werde noch auf y, = 0 normiert. 
Das für die Differentialgleichung (4) zu lösende Randwertproblem lautet 


G=0 füry=0, also ©,(9,0) =0 (y = 0 ist Düsenachse), ) 
9=e (y)fürp=0, also @,(0, 9) = e(y) (gegebener Einlauf bip=0), | (6) 
= 9,(p) für y=y, also @,(9,%Y,) =9,(P) (y = y, ist Düsenkontur). | 


Man setze nun für das Folgende @(g) als von einerlei Vorzeichen voraus?). Dann kann man in (4) 
neue Veränderliche ®, %, 2 einführen durch 


RN 
= , -. NE MB, Many)... ne. (7). 
Dabei soll noch a höchstens integrierbar unendlich werden für 9 =0. Dadurch geht die 
Differentialgleichung (4) in die Wärmeleitungsgleichung 
DI FI een (8) 
für 2(®8, Y) über mit den Rand- bzw. Anfangsbedingungen } 
2y(8,0)=0; 0, P)=EP; Bd, P)=KRP)...:...- (9). 


Hierbei entspreche %, dem Wert y,, R(®) sei das von p auf D umgeschriebene ©,(p) der Vor- 
gabe in (6) am Rande der Düse, und es sei E(Y) = e(y) die Vorgabe am Einlauf. 

Da mit 2(8, Y) auch Ay(D, Y) Lösung von (8) ist, kann statt des obigen Randwert- 
problems für 2(®,Y) auch das folgende Randwertproblem für ©(®, Y) = Qy(®, P) formuliert 
werden: Eine Lösung von 


3 nu. (10) 
| 0 zu finden, daß die Randbedingungen 
« 9(8,0) — 0 für 20 
€ 90,9, =E(P) für OSY<SP, (mit E(0) = 0) MEER 
.9(8, W,) = R(®) für © 20 (mit R(0) = E(%,)) 


* erfüllt werden. In der Theorie der Wärmeleitung entspricht dies dem Problem, die Temperatur- 
_ verteilung in einem endlichen, thermodynamisch homogenen Stab als Funktion des Ortes und 

der Zeit zu ermitteln, wenn zur Zeit # = 0 die örtliche Verteilung und wenn an den beiden Stab- 
enden die Zeitabhängigkeit der Temperatur gegeben sind. 


. ) Hieran liegt es, daß nur der unsymmetrische Schalldurchgang (Meyerscher Typ) behandelt werden 
oht aber auch der symmetrische (Taylorscher Typ). 


Au 


1 
ir. 
ng 


u de A a a ni Pal und 1.0 ee 1 A EEE EL ZN En -. 
“ e - Pa N ' - 
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3. Mansetze R(®) für D® >0 zweimal differenzierbar mit absolut integrierbarem R’’(®) 
und E(Y) in eine konvergente Fourierreihe entwickelbar voraus. Sei 


> 14 | 
EP) = 3 6,sin ng 0 (1) 
= j 


diese Fouriersche Entwicklung. Dann wird das durch (2.10) und (2.11) gegebene Problem für 
8>0,0sYSsY, durch die Funktion 


) 
v Ss sinnz -u= > 

9 {i et zei 5) “RN n 1 En 4 7 

ao IE u Ro+|e zoa 3 
{}} 

wo 
+2 sinn de e 
gelöst. 


Mit den getroffenen Voraussetzungen lassen sich die Konvergenzbeweise für die auftreten- 
den Reihen unter Zuhilfenahme des Abelschen Kriteriums für O <Y<Y, undalle © 20 
leicht führen. Die einmalgliedweise nach ® und zweimal gliedweise nach Y differenzierten Reihen 
konvergieren gleichmäßig für 0 sY<WY, und alle © Zö bei jedem 0 <Y%,<Y, und jedem 
5 >0,sodaßfürO <Y<Y, und ® > 0 die partiellen Ableitungen O8 und Oyny durch glied- 
weises Differenzieren gewonnen werden können. Damit läßt sich die Tatsache, daß ©(d, %) 
gemäß (2) Lösung von (2.10) mit den Anfangs- und Randbedingungen (2.11) ist, leicht nach- 
weisen. Eine eventuelle Ausnahme macht nur der Punkt © =0,Y=W.. An dieser Stelle kann 
die Einlaufsbedingung nur dann auch noch erfüllt werden, wenn R(0) =0, d.h. Parallelität 
der Düsenkontur zur Düsenachse am Düseneinlauf gefordert wird. 

Geht man in (2) von ®, Y wieder zu 9, y zurück, so erhält man dadurch die entsprechende 


Darstellung für den Strömungswinkel © im Felde der , y 


Y 22 Pn 2 nr? 
.. Sam ae = MER 
(+1) y? L 2 r 
99,1%) = + 9,@)+- > (1 — a, rs | O4 Jet Ride 
Yı 7 n=l N ‘ (3) 
0 
Sera 
oo @t) 
I sin nnte ) 
Nn=1 Yı 


und hieraus schließlich auch vermöge (2.1) und (2.8) 


2 
W=W(,y) = en Bi Ba gend 


d.h. schließlich 


W(o,y) = — 
liege 
wo ©, durch (gliedweise) partielle Differentiation aus (3) zu gewinnen ist. 
. 4. Schon die einfachste Wahl, nämlich 


G(p) = const = AUS Free ee 2 oe (1) 


führt zu einem praktischen Beispiel für die allgemeinen Entwicklungen, das alles physikalisch 
Wesentliche zeigt. Wegen @(p) >0 wird die Düse vom Medium .mit wachsender Geschwindig- 
keit- durchströmt. Es ist 


— =G(9)=R, also W(p) = — P*). 
Damit folgt aus (2.5) bei y, = 0 


9,(p) = (#+1)y1 At(@ — Pp*). 
Nach (3.3) ergibt sich daher 


[3 & y 
8 nn sinnn— - 
O(9,y) = (+ 1)Ay(p — pr) +2 (x +12 46.32 (— 1)" Pr 

U pen], n3 

nn? ne? 
oo 2 ee a 
+3 dusinnnt re man 
ne, Yı 


E“ 


+ e .* . pi 
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mit (— 1)" («+ 1)y? 22 
da er rl 5; (2), 


oder durch Summation der ersten der beiden Reihen rechts 
‚3 gr 
N ee 
y 


© SGrDRE 
pp) = (+ 1A EEE eye) HEdsinnnd. TR 0) 
n=1 N 


Aus dieser Verteilung der Strömungswinkel © in der 9, y-Ebene erhält man gemäß (3.4) die 
Geschwindigkeitsverteilung 


x Y * 22 
W (o,y) = Bi; een] NREEIELr 3 nd, cos une 2 4 wrDR N "a 


ale 
6 Yı ae1 
in der 9, y-Ebene. 
Alle durch (3) und (4) gegebenen Paare von Funktionen W(o, y), © (9, y) liefern Düsen- 


strömungen mit konstantem Gradienten @(p) = der Stromfadengeschwindigkeit im oben 


d 
dp 
definierten Sinne. Ihre praktische Brauchbarkeit ist wegen des näherungsweisen Ersatzes des 
Systems (1.7) durch das System (1.11) allerdings auf die Umgebung des kritischen Strömungs- 
zustandes beschränkt. 

Unter den in (3) und (4) gewonnenen Lösungen zeichnen sich diejenigen besonders aus, 
in denen d„ = Ofürallen =1,2,3...ist. Dies macht es erforderlich, daß die Fourier-Koeffi- 
zienten c„ der den Einlauf charakterisierenden Funktion e(y) gemäß (2) durch 


Y 3 =, “+Dy, | 
—=2— A Era th Sue un 28T, 5 
le) + () 
bestimmt werden. Ein solcher Ansatz ist aber erlaubt, da für die durch (5) definierten c, die 


unendliche Reihe 3 c} konvergiert, die c„, also nach dem Rieß-Fischerschen Satz tatsächlich 
n=1 

die Fourier-Koeffizienten einer quadratisch integrierbaren (im vorliegenden Fallsogar im Inter- 

vallO<y<y, stetigen) Funktion e(y) sind. Damit hat man dann die einfache und elementare 

Polynomiallösung 


1 
W(e,y) = 29 — p*)+ 


May? — y)) | 
1 
9p,y) = @+ A pp a ae Zn 


gewonnen. Sie zeigt schon alles Charakteristische einer Düsenströmung in der Nähe des Schall- 
durchgangs und soll daher im folgenden noch etwas näher diskutiert werden. 


Für die kritische Isotache W = 0 ist nach (6) 


1 
AypP—Y). 

Ihr Schnittpunkt O mit der Kanalachse sei 9 = 9,, so daß 

a 


0 = (pP — gr) + 


=) Ay: 
gilt. Für das Folgende werde O als Koordinatenanfangspunkt gewählt und noch | 
Berenwey . 3. 1 Dei... (7) 
gesetzt. Dann gehen die Gleichungen (6) in 
kin +1 _ 
W=# 4% | 
EB.) RT 
9 = (x +1)99 + en 78 


über. en Ursprung O der 9, Y-Ebene entspreche der Ursprung © = y = 0 der x, y-Ebene. 
Erueßlieh führe man noch die Ähnlichkeitstransformation 


or 1 (9) 


4 aus und ordne die 9, 9-Ebene der &, Y-Ebene zu. 


a 


ve 


} } = F 5 | u 
Kt | 108 Behrbohm, a d es ul unsymmetrischen Schalld 


Dr 
Be Ron 
R Die Isoklinen © = const sind im allgeme Fall Kurven dritter ng der 9, ig: bene. 
4 3 Nur die zu9=0 gehörige Isokline zerfällt einmal in die Y-Achse, das andere Mal in die durch, 
> Re - N 4) 
E n = yp? . . RE . Pa: Swirl! Dr = 10) 
gegebene Parabel. Dieser nichtlineare Zweig der ni 0 =0 ist der DE ische Ort der. e 
achsennächsten Punkte der einzelnen Stromlinien. Er ist in der &, 2 -Ebene nach (1.4) durch Ben, 
»HIsoe "RE 
N‘ = — 5 ydy, dy= E; 
® c* o*c* Ä 
gegeben. Hierin sind ”s und Tr als die gemäß NR und (10) durch 
W<- & x nn 1 m 
Bu it ; bestimmten Funktionen von 9 aufzufassen, so NL man damit dann die Darstellung = 
BIS de | 4 Pr wi; ; 
= % EINER 4 a ? 
In = 773 IEre. : 2 RE ee 5 nee (11) 


für die Koordinaten 7,, Y„ der achsennächsten RN, mit Hilfe des Parameters 
erhält. 


Für eine Stromliniey = const =Y, An 


‘ 


Be 


2 w RE AA 
da cos Od9, dy=—sin Odp, 


wo W und © die nach (8) durch 


Ze) TE BF DIE 2 er 
wm oa... 
definierten linearen Funktionen von 9 sind. Damit wird | 3 
ur. o N o N u z P> 
E; RT or = Baar. OE, = RS”, 
N. z = my) + | „cos Odp, Y=YmYı) + IF ER nl: 
Satz Be | 


die ee = (9,9 = 36) der Stromlinie 9 v= Y1, wo in). Ym (91) nach.) 
© nach (12) und — — über (1.10) nach (12) zu bestimmen sind. 


ö ' Die Gesch einige Ver auf der Düsenachse (Y= 9 ist nach ( (12) durch WM. =. 
‚also implizit durch 


15- a5 -Zaw-® 
gegeben, woraus (bei x = 1,4) 
wje ı(®) / 
SR oe \c* -z ne ae a + 2arc tg dr 


Bl l/wie/ 9 ei 1,44 (85 
Foyer “ Bun 
‚folgt. he: ; 


Schließlich interessieren noch die che W=- con 
die Parabeln 


gegeben, so daß für sie 


A 
“ 
h, 


[7] N m N % h, « 
z 4 ER EN > 
€ - a = 4 n 
T h € ne EN 
' r , 
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; gilt, Nach (d 4) Abd ( (16) ist EL 2 
N c* 
di= -efer 1) y cos + "sin 0 dy, dy -" (+1) ysin 04.005 9)ay (17). 


t (15) ist © nach (8) is W = const die durch 
0-4 


definierte kubische Funktion von %. Hiermit und mit dem aus (14) zu W =const bzw. dem 


gemäß (1.10) zugehörigen er bestimmten © = z, erhält man die Isotache W = const aus (17) 
durch Quadratur in der Parameterdarstellung 


Le, GE r * 
D ee &(Y) =. | (mt FO +2 sine) an 


N c* f ” o* a 
y=y(%) = (41950 +8 cos 0)d7 
ö 


Schließlich noch ein Wort über den sekundlichen Massenfluß durch die durchy =yYy, 
| Done 
festgelegte Düse. ‘Der Definition der Stromfunktion zufolge ist dieser durch 0 = - % Y, ge- 


; * 
geben. Gemäß der Stromfadentheorie wird er näherungsweise durch Q; = = - Yın gegeben. 
f Yı 


0o*c* 


"Wie der Vergleich mit (11) zeigt, ist dieser Näherungswert etwas zu groß, dain ym = br dy 
o* cr 
der Integrand = ——, die reziproke Stromdichte, größer als list. Die Korrektur, Sek auf 


den Näherungswert, ist durch 


u. NE A ER A TEN 19 
for = (19) 
0r.0* Br 
dy 
ol ER de 
gegeben. 
Für die Krümmung k der Stromlinien gilt allgemein 
a, _ 0 de uae 


a > 4 
und daher für y=%, nach (12) 
ö en \ : Sau}, 
k=@+H) Rn 


Be e: hdsten Stelle: ist W.= ee, v woraus sich — —, über (1.10) und damit % = / an der 


5 
5 
„ 
us . 
te 
ir nn ,% 
x . > 


9 


Stelle beim, Für den tuning na Rn = m an der engsten Stelle gilt daher 
3 AR. Br 
I ar 


m 
s he Verhältnis von Halver Düsenhöhe zu Kr ümmungsradius an 


ee 
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Bild 2 und Tafel 2 orientiert über den Verlauf von « als Funktion von y,, Bild 3 über die 


Geschwindigkeitsverteilung auf der Düsenachse. 


a ee‘; 


N 
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I 
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Tafel2. 2 
06 
q = = 
PR a Yı Yı a 
Pr 0) 0 0 0 
0,1 0,2006 0,04 0,0038 
Br 0,2 0,2792 0,08 0,0154 
0,3 0,3359 0,12 0,0350 
0,4 0,3818 0,16 0,0627 
z 0,5 0,4203 0,20 0,0992 
0,24 0,1456 
0 0,28 0,2010 | 
0,32 - 0,2683 j 
TEE 00 | 040 
0,44 0,5551 [ 
In Bild 4 wird ein Bild der Strömung gegeben. Gezeichnet 


Ym 
Biid2. Krümmungsverhältnis a = —— als : a 
wurden die zua =0,1; a = 0,3 und a = 0,5 gehörigen Strom- 


BANK MOR ae bei der Jinien. Diese drei Werte wurden gewählt, um in den beiden 
— halbe Düsenöffnung an der engsten ersten Werten den Anschluß zu gewinnen an die von Oswa- 
Pe - titsch und Rothstein[4] berechneten Hyperbeldüsen und um 

Am nen an ‚der OngBien” nn ELEHSNVert (@a= 0,5) wesentlich über die Oswatitsch- 
Rothsteinschen Rechnungen hinauszukommen. Während näm- 

lich das dort entwickelte Näherungsverfahren an kleine Werte von « gebunden ist (dafür aber 
im ganzen w-Bereich gilt), ist eine derartige Beschränkung für die hier vorgetragene Lösung 
nicht notwendig (sie gilt dafür aber auch nur in hinreichend enger Umgebung der Schall- 


Ym 


13 


e* 


z 07 Er 


Bild 4. Stromlinien und Isotachenbild der Strömung (4.8) 
—— — — Isokline © = 0 


Bild 3. Geschwindigkeitsverteilung auf der Düsenachse 
„der Strömung (4.8). = =2:a; Y = A:y mit A2 = G(p) 


geschwindigkeit). Außer den genannten drei Stromlinien wurden noch die Isotachen w/c* = 0,8; 
0,9; 1,0; 1,1; 1,2 und 1,3 berechnet und eingezeichnet, überdies die Isokline © = 0, der geome- 
trische Ort der achsennächsten Punkte der einzelnen Stromlinien. i 
Os rg. als 


Bild 5 und Tafel 3 gibt schließlich die Mengendurchflußkorrektur Q 
s 


Funktion von a wieder. Die Oswatitsch-Rothsteinsche Näherung 


0 —Q SS a? 


ee 


Qs 90 


wurde gestrichelt mit eingezeichnet. Es zeigt sich gute Übereinstimmung. 
5. Die aus (4,6) gewonnene Lösung (4.8) des Systems 


20 oW 20 oW 
ae el — =(2+1) W — + 

op 9% 9y Ip 

zeichnet sich nicht nur durch ihren elementaren Charakter, sondern auch durch eine weitere g 


wesentliche Eigenschaft aus, die bisher noch nicht besprochen wurde. Für ihren Gradienten W, 
bzw. W; gilt nämlich A ar i 


W, = 4: bzw. W,=1, 


an a ai Dr ea PR En 1 A a A Fra Kr ES TE TE A SE U FA I u 
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d. h. hier ist nicht nur der aus der Stromfadentheorie entnommene Näherungswert des Gradienten, 
d.h. W,= = 4? = const, sondern sogar exakt der Feldgradient W, = 42 = const. Aus diesem 


Grunde kann die gefundene Lö- Tafel an 
sung (4.6) in dem Sinne als sehr 
gut angesehen werden, daß we- a (Qs—Q)/Qs 

mipstenssein Teil der geirallenen Sn | Er 
Näherungsannahmen, nämlich 0 0 

der auf den Ersatz von W, en 

durch ga in der Differential- 0.0992 —0.00021 -0003 

gleichung (2.2) bezügliche, hier- 0,1456 —,000051 

bei ausgeschaltet wird, so daß N BR at 

nur noch die Annahme des Er- 0.3482 _0.00301 Bi 

satzes von K(W) durch seine 0,4430 —0,00481 

Tangente im Schallpunkt W = 0 0,5551 —0,00744 t 9006 


sich auswirken kann. 
. . . . . £7 
Mit dieser Erkenntnis erhebt sich daher sofort die Frage | 407 
nach weiteren Lösungen, für die W, = W, im ganzen Felde 
gilt. Zunächst werde die spezielle Frage nach weiteren Lö- Bild 5. Mengendurchflußkorrektur als Funk- 
3 7 tion des Krümmungsverhältnisses an der eng- 
sungen mit G(p) = W, = W, = 4? = const behandelt. “sten Stelle für die Strömung (4.8) - 
Aus W == 22 folst W 2 also W. En Q= exakter Mengendurchfluß 
DE 16) d ah fe) zolg .s rip» ie 16) Ba a (Y) Qg= Stromfadennäherung der Durchflußmenge. 
N ee =f'(p)p +g(y). Daher ist ©, = f’(Y)P ___ 1. Nanerung nach Oswatitsch-Rothstein. 
+9g(y), was gemäß der zweiten Gleichung (1.11) zu der 


Forderun u { 
i Pt) = +1) (pt SW) 
führt. Der Koeffizientenvergleich in p führt dann zu den beiden Forderungen 


FW) =(k+1)A, gp) = (#+1)2f(y) 
ar 


woraus 


|’ = May? +a,y+a, 


\sn- «+ Fory+ rat at lRar+a, 


folgt. Bleibt man bei der ne in Ile auf eine Achse symmetrischer Düsen, so ist 
a, = 0 zu setzen, da nur dann W eine in y gerade, © eine in y ungerade Funktion ist. Damit 


hat man dann aber, a, = u 2 — 2 0*, a, —=0 gesetzt, gerade die Lösung (4.6) gefunden, 


was besagt, daß die Lösung a die Ei ihrer Art für W, = const ist und daß sie daher 
durch diese Eigenschaft charakterisiert werden kann. Daß nämlich die Annahmen a, = a, = 0 
keine Beschränkung der Allgemeinheit der Betrachtungen bewirken, sieht man leicht so ein: 


Ads = 0 ist sofort durch Addition einer geeigneten Konstanten zu © zu erreichen, was gemäß 
2 


(1L.11)ohne weiteres gestattet ist; «, — 0 kann man durch Übergangvon 9, yzu 9, =P— DIPzuT ; 

y,=vyY+ is erreichen, was aber nur auf eine andere Zählung von @ und y hinausläuft. 

Im allgemeineren Fall, daß an der Annahme W, —= const nicht mehr festgehalten wird, 

folgt aus W,=@(p) =H"(p) sofort W=H'(p)+f(y), also 0©,=W,=f’(y) und daher 

O=f’(y)p+g(y), woraus sich weiter ©, =f’(y)p+g’(y) ergibt. Die zweite Gleichung 
(1.11) SR damit zur Forderung 

Wet) =RtHYD(HOHSW) EM) 2: (1). 


Hieraus folgt durch zweimalige Ableitung nach 9 
[H’(p) +] HP (9) +3 4" (p) HE) = 0, 
woraus sofort f(y) = & = const und damit nach (1) auch we — ß = const folgt. Das heißt aber 
Bel ENTE Se RO ER (2). 


Für $ = Oist Hann entweder Z"(p) = 0, also le —O, oder Er p) = — a, also wieder G(p) = (), 
welcher Fall zum Potentialwirbel führt. 

Ist aber $ #0, so erhält man aus (2) 

d 


Er (H’2(9)+2&H’(p)) = 1 
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und daher bei beliebigem konstantem y 


Wi =—at) tr +Y 


eh 


so daß man hiermit auf die Quellströmung 
Ww=+) 
geführt wird. 


5 
—ptr, O=BuHb 
Leitet man (1) je einmalnach p und y ab, so wird man durch die entsprechende Diskussion 
zwangsweise auf H"(p) =4? und damit auf die Lösung (4.6) geführt. Die einzigen Lösungen 
des Näherungssystems (1.11), für die der Geschwindigkeitsgradient W, nur von p abhängt 
und dann natürlich dem Stromfadengradienten W, gleich ist, sind die Potentialwirbelströmung, 
die Quellströmung und die durch (4.6) gegebene Strömung. 


| 
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Zur hydrodynamischen Schmiertheorie des Zapfenlagers 
Von ©. Weber in Schlewecke 


„Strebt nach der besten Lösung“ 


Ausgehend von der Grundgleichung und vom Variationsproblem des endlichen Zapfenlagers wird die 
optimale Näherung der Druckfunktion gefunden für den Fall, daß die Näherung als Produkt einer Funktion 
von x (in Laufrichtung) und einer Funktion von z (in Achsrichtung des Zapfens) gewählt wird. Für die 
Funktion von z finden wir eine geschlossene Lösung, für die Funktion von x eine Differentialgleichung. Für 
diese wird an Hand eines Beispiels eine Lösungsmethode dargelegt. Für das Beispiel wird die Lagerkraft 
und ihre Richtung berechnet. 4 


Based on the fundamental equation and the variation problem of the finite pivot bearing, the best approxi- 
mation of the pressure function is found out in the case that Ihe approwimation is composed of a function 
of & (that is the direction of motion) and a function of z (the awis of the pivot). The latier is given in a closed 
form, and for the function of x,adifferential equation is stated that is solved. by a special method in an example, 
” " The example is extended to Ihe computation of the quantity and the direction of Ihe force in the bearing. 


„En se basant sur l’öquation fondamentale et sur le probleme de variation du palier fini d’un tourillon 
on trouve la meilleure approximation pour la fonction de pression au cas que cette fonction soit le produit 
d’une fonction de x (direction de mouvement) et d’une fonction de z (direction de l’axe du tourillon). Celle-ci 
est assignable sous forme finie, tandis-que pour la fonction de x une Equation differentielle est avanc£e, dont 
la solution se fait dans un example par une methode speciale. Pour cet example quantite et direction de la 
force sur le palier sont determines. A 


Mexons u3 OCHOBHOrO ypaBHenusa u u3 BAPHAILHOHHOH 3aayH IA TOMMIHIIHURA koHeyHohi ; 
f . VIMHBI ABTOP HAXONUT OINTHMANIBHOe IIPHÖNMKEHHOe PenıeHue AA byHKUHM naBleHHaB 
TOM Cyyae, ECM 9TO IPNÖNMKEeHNe MPENCTABIAeTCH B BuNE IpomsBeneHnn DyHEumH OT 
© (B HAINpaBıIeHnn BISKeHNT) Ha PyHEIMO OT 2 (p Hanpasıenun ‚ocu). Ilocıenuas gaerca 
B 3&MKHYTOM BuNe, a IA PYHRIUM OT x HAXONuTCA MuhepeHumanbHoe ypapHuenne, MeTOoA — 
PpemeHnsd KOTOPOrO NOACHACTCH Ha Ipumepe. lust Toro npumepa oMpepeuserch ONOopHaA 
CHA IIO BEIMYHHe NM HALPABIEHNID. 42 Na Br 


d 
’ 


A. Grundgleichungen 
Für das Zapfenlager nach Bild 1 erhält man nach der Reynolc 


‚für den Druck p die bekannte Differentialgleichung 
/ d 3 d * UF) 3 2 . ru 

wand, 

9x\n dal 9z\n oJ) ° 
Hierin ist 4, die Umfangsgeschwindigkeit des Zapfens, h die 
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Schmiermittels, & die Koordinate des dünnen Schmierfilmes in Bewegungsrichtung und z die 
Koordinate senkrecht dazu, also in Achsrichtung des Zapfens. Die Randbedingungen lauten: 


B 
»=0 für z2=- — (Seitenränder) 


2 
= 0. für &=%, (Einlaufrand des Schmiermittels) 
0 
» = 0 Zi = 0  (Auslaufrand des Schmiermittels). 


Gl. (1) ist die Euler sche RT ne folgenden Varjationsproblems: 


3 9 3/9 . 


Die Zähigkeit 7 sei & ie ee, Unter- 
suchung konstant angenommen. 

An Stelle von & führen wir die Winkel- 
koordinate @ ein: 


I, 


p. 
Hierin ist R der Lager- und r der Zapfenhalb- 
messer. 


In gleicher Weise führen wir auch- für z 
eine neue Variable ein: 


R-+r 
2 


R-+r Bild 1 
2 = ——vy. 


2 


Die Randwerte von z sind SeBoN 2 ‚ die entsprechenden Randwerte von y sind + ß, also 


Die Spaltstärke wird, wenn wir p von der engsten Stelle an messen, 


ee 7 


2 


..R—r ’ N i i HERR Pe e i 
Das halbe Lagerspiel SIERT setzen wir gleich s; eist die Exzentrizität des Zapfens, „bezeichnen 


wir mit y, dann wird . 
h=s(l—xcosp). 


Nun führen wir an Stelle von % und p folgende dimensionslose Größen ein: 


h*= L = (1—yxcosp) 
nd nn ENDE, 2 REES PR, iR. (9). 
a 
Nu 9 
Dann lautet die Differentialgleichung 
d en) 02p* . 0h* 
De re ES Bat —— el een rl 2 ( 
| | zu TI . 


mit den Randbedingungen 


p=) fü 9=M 
3 ar. 
ME 0 für dp N) 
Das dazugehörige Variationsproblem wird: 
9p* 9p* 9p* - 
h*3 er -) + An )- 2h*—-dody=Min ....... OR 
eu Oy op a 
8 


{ un 
RR 3 iR e i „6 Ks 
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B. Optimale Näherungslösung pt =f()e)) 
Wir suchen eine Näherungslösung für p* von der Form M- 
u, Pr=N)IP) - nern. 
Ba, Die Randbedingungen lauten hierfür 
| g(+P) =0 
Be, go) =0 a afio,) 
BEFNPL 
72 — 0 PERL m 0 
Y J(pı) do 
B: Die Anlaufstelle ist für diesen Näherungsansatz zu einer geraden Linie  =9@, geworden. r 
e Für p*=fg wollen wir die optimale Lösung finden, für die, 2 
Bir die linke Seite der Variationsgleichung (ö) zum Minimum wird. 7 
Br (Das hierbei angewandte Verfahren läßt sich mit Erfolg auch bei anderen gleichartigen Pro- 3 
2 blemen anwenden.) df(e) A 
Ro. Hierzu führen wir in Gl. (5) für p* den Ausdruck nach Gl. (6) ein, setzen er. = fund RS 
: dgy) dp | 
—— = g’ und erhalten: 
day ae | i | | ; 
T A fg hr — 2 Ari gldodp—= Min. 2... (0% ee 
—ß 9 F 
Wir wählen nun ersteine beliebige Funktion f (9), die die Randbedingungen 4 


für 9 = 9, und p,erfüllt, und untersuchen, welche Bedingungen die Funktion g (y) erfüllen muß. 
Ist f(@) gewählt, so kann man in Gl.(7) die Integration nach p durchführen. Die hierbei 
auftretenden bestimmten Integrale bezeichnen wir wie folgt: 


Te 
RE 
P 


9 m ii 7 
a, — [h*f’?dp Be a ae de he 20 Ge 


Po 

Pı 
As — [hr f2do 
b Po 
Dann erhalten wir aus Gl. (7): 
+ß j SEIT ser \ 
J, [y®+ag?—agldy=Min.... 2.2.2.2... 
Hieraus folgt für g(p) die Differentialgleichung | 


EL » FRE oe ae ar a (10). A 1% 

(da 4, \ EEE 

Die allgemeine Lösung dieser Gleichung lautet bei Beachtung der Randbedingungen g(+ Pf) = 0: 

an {A EEE > 

mh Eier N 
ie | en Cojapl EEE bl Ei RT une 
‚ mit a2 — = und ig —— ®n I h “ ii; vn, Y 

Ay dr) Pr 


Um: die optimale Lösung für p*(p,y) = f(g) y(y) zu erhalten, muß für jede beliebig ge 


a» Dun 


N 5, 


RE 
A F 


B . LI er “ 
Jg" (V) Au Eu 


!) Ich verweise hier auf folgende frühere Näherungsansätze: 
des Berl. Bezirks-Vereins dtsch. Ingenieure, Nr. 5, Mai ana Tone 16 
unendlich langen Lagers und für g (y) eine Parabel. — G. Vogelp 
3. Heft 1943, $. 192—212, ebenfalls für f(p) die Funktion des unendlich 
_ die optimale Lösung mit Hilfe einer Reihenentwicklung. Diese Jäl t si 
gezeigt wird, in geschlossener Form ausdrücken. a , 


’ » Pr in 


T ® DR x 1% 
SUR a Fr BE 


%. Anden. Math. Mech. 
Bd. 30 Nr.4 9 


Ne. u | 
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-Nunis it folgendes zu beachten: Die Näherungslunktion 2»* (o, y) ist das Produkt 
der beiden Funktionen f(@) und g(y). Bei der Zerlegung der Funktion 9* in diese zwei Faktoren 
können wir die eine Funktion mit einer Konstanten multiplizieren und die andere Funktion 
durch dieselbe Konstante teilen. Das Produkt p* (p, y) ändert sich hierbei nicht. Multiplizieren 
wir 9(y) mit 2 
1 

ten Integrale des Gleichungssystems KE denselben Wert. Wir können folglich unbescha- 

det der Allgemeinheit annehmen, daß g(y) so gewählt ist, daß b, = b, wird. 
Führen wir b, = b, und b, in Gl. (7) ein, so erhalten wir 


und bezeichnen wir — m wieder mit g(y), so geben die zwei ersten bestimm- 


$ı . 
[rer +2 ws anf eu 2.2013). 
i 
P 
Hieraus folgt die Eulersche ee Br F (9): 
LE . 1x | 
Gh do TE !- ee... da, 


Eine Lösung in geschlossener Form ie, wir für diese Gleichung nicht angeben. Im 
nächsten Abschnitte wird an Hand eines Beispieles ein Verfahren gezeigt, wie eine Lösung für 
f($) gefunden werden kann. Um eine solche Lösung zu erhalten, sind folgende Werte anzunen- 
men: 


b, 
LT. DE Wert = wir bezeichnen ihn mit I ) s 
a 


. Der Wert für 4; damit ist A* =1— x cosg als Funktion von @ festgelegt. 


3. Der Wert 9,, für den f(p,) und Zn gleich Null werden. 


Mit der Annahme dieser Werte liegt der Verlauf von /(%) fest, wie das Verfahren zeigen 
wird. Die Funktion f(p) wird für verschiedene fortlaufende Werte o rechnerisch bestimmt, bis 
sie wieder den Wert Null annimmt; dieses gibt den Einlaufrand @ = 9.. 

{ . a bs 
Wir wollen annehmen, daß wir für die drei Werte (%) ‚x und o,den Verlauf von f(p) ge- 
A: 
funden haben; nun müssen wir die dazugehörige Funktion g9(%) finden. 
Hierzu bestimmen wir nach Gl. (8) die Werte a,, d; und a,. Zwischen diesen Werten und- 


b j I EraE > 
dem angenommenen Werte () besteht eine Beziehung, die im weiteren von Bedeutung sein 


k Ya. 

wird. Um diese Beziehung aufzustellen, gehen wir von Gl]. (14) aus. Durch die gefundene Lösung 
/(g) wird diese Gleichung zur Identität. Wir multiplizieren die Gleichung mit f(g) und inte- 
grieren sie über @ von 9, bis 9,. Dann erhalten wir 


Pı 
*3 er EN () 3 f2 dh* Ai 
J B (1 de f hr Bi PB — 75 dp =0. 


= Das erste und dritte Teilintegral wird partiell integriert. Die Randintegrale werden gleich Null, 
da f(p,) und f(g,) gleich Null sind. Die nachbleibenden Integrale sind die von uns eingeführten 
Werte a,, @, und a,. Wir erhalten hiermit 


b r 
+), = ee ea ir (15), 


Nachdem a,, @, und a, bestimmt sind, kann die Funktion g(y) nach Gl. (11) gefunden 
werden. ä 

Ei Bestimmen wir zum Schlusse aus g (y) nach Gl. (12) die Werte d,, 5b, und b,, so müssen 

folgende Bedingungen erfüllt sein: 


1% IR Pie 


ri: br (b 
€ Ar. 5% ar: ;). ö 


In Wirklichkeit handelt es sich aber nur um eine Bedingung. Um dieses nachzuweisen, 
ehmen wir die Bestimmungsgleichung für g(y), Gl. (10). Diese ist durch die Funktion g(y) 
t, so daß sie zur Identität wird. Wir multiplizieren die Gleichung mit g (y) und integrieren 
Be vi von —ßbis +ß: +5 


n a a | £ 
f 0 g+lglay =0. 
> I2 Pu h 
‘ g* 


u ARTE 
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Durch partielle Integration wird diese Gleichung umgeformt. Führen wir nach Gl. (12) die Werte 


” ” ” . 23 
by. b, und b, ein, so erhalten wir nach Multiplikation mit 2: 
j 1 


b, 
Pe el Sand BE 2 5.2 An a (16). 
1 


Ist die erste Bedingung 5, = b, erfüllt, so folgt aus dem Vergleich der Gl. (15) mit Gl. (16), 
daß = = ) und somit auch die zweite Bedingung erfüllt ist. 
a 
Wir müssen also untersuchen, unter welchen Bedingungen b, = b, wird. 
Wir nehmen hierzu die Lösung für g(y) nach Gl. (11) 


nn a a 41 , % 
mc) Te, Te 0 Be und en; 
Hiermit finden wir b, und 5b, nach 
aß} Gl. (12) ß 
4 
2 {- +B . N 
EN LIE 2 [eh-g 2] 
Zu &oj 
4 + ——; — au ee) mr | 2 —_ 
35 a | a BER el — Ba c- [& 55 & ß] 
; IR an +B e 
25 \ 201 &Yy 
eile f Bi la 
ABER GERIET IT IE Enmie 2 
a 1 
7 — (2 zu A 3Lg aß + Gera]: 
i ST te 
1-1 1 = Ei Aus der Be b,=b, folgt hieraus 
GNS RTL TE EN 225.30 5 RE E31 BR S 2 (aß —Tgap) MER N 
iR RERE aß 
Bild 2 2aß —3Tgaß + paB 


Da € bestimmt ist, so kann nach Gl. (17) der Wert für «&ß, und, da auch & bekannt ist, der Wert 
für ö gefunden werden. Zur Erleichterung sind für einige Werte «ß die dazugehörigen Werte C 
berechnet und iin Tabelle 1 angegeben. Weiter ist der Verlauf von «ß als Funktion von Cim Bild 2 
dargestellt. 


Tabelle 1. 
RR. |, 1,07" 11,07 RE a 
CI =11,01: 3,83. [72,140] 1,66 2 1 a 1 DES DEE ee 


. 


C. Numerische Näherungslösung der Differentialgleichung (14) 


> die Funktion f(p) zu ermitteln. 


Aus dieser Gleichung ist für einen gewählten Wert 
\ ; 1 
Außerdem sind noch Werte für y= — und für den Endpunkt 9 =9, zu wählen. Wir zeigen 


hierfür an Hand eines Beispieles ein allgemeines numerisches Verfahren. Für das Beispiel 
grählen wir folgende Werte 


Oo 


Del ug =08, 


Ri a 
Wir beginnen mit den Anfangswerten /(9,) =0 und u @.) =0. Dann berechnen wir R 


Busreise hieraus die Werte f(p) und Te für p = 20°, 150 Ne Die Intervalle sind folglich 

gleic a | 
Ag = 5° = 0,0873. 

Die Zahlenwerte, die bei der Berechnung auftreten, sind in der Tabelle 2 eingetragen. Die 

erste Spalte gibt die Werte 9, die zweite Spalle die a Werte für cos p und Sin p. 


Ei dh* 
Nun müßten wirnoch A*=1 2 c05p und — = # sin 1p eintragen. Um die Rechnung zu ver- | 


dp 


DE 28 1. a an ie rt 
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. i 1 1 * 
einfachen, nehmen wir statt dessen - = —- — cos =2 — cospund B di =sin 9. Dieletzten 
X 4P 
Werte sind schon in der dritten Spalte enthalten, die Werte — h* sind in der vierten Spalte, in der 


: X 
fünften Spalte die dritten Potenzen = h -) eingetragen. 


: j & d d 
Nun finden wir den Mittelwert von ar il für das erste PuSV allvong= 20° bisp = 25°, 
d 


Um wieder Rechenarbeit zu sparen, berechnen wir statt dessen — 212 Ras A er 
een GO Summe: 4 \do dy/mitteı. Dafür setzen 


fl: —(d = 1 -( = 1 AH =  (£ dh* 
h*3 +— we de 2 2 7x3 h*3 Fe N 
dp dl g=25° dp dp 9-20 ur ir dp /g=25° 2% v4 ‘+ dp /g=2* 


Hierbei ist f (25°) Be Null, f (20°) hingegen noch unbekannt; wir setzen es in erster Nähe- 
rung ebenfalls gleich Null und verbessern die Berechnung nachträglich. Für das Intervall von 
p = 20° bis @ = 25° wird somit 


(ur x 4 zZ) = sin 250 - sin 20° = 0,7646. 
dp dr mittel 


e ec e Na: 

Wir multiplizieren diesen Wert mit 32 do =2 Ag = 0,1746 und erhalten den Zuwachs von 
= h*3 an ‚ also 
2.2.49 


v 


Re A (ms =) = 0,1838 
dg 


m 
Da Ind für 9 = 25° gleich Null ist, so wird 
ER dpa | 


4) % 2 
(Gi Tl 0 1335- 


*\3 ; 
Nach Tabelle 2, Spalte fünf ist 1) —= 1,192, also wird (2) = — 0,1335 : 
1,192 = — 0,1120. X g=2° dg/g=20° 
df 


Der doppelte Mittelwert von z wird für dieses Intervall 
[0 = 


= in 
a — 0,0 — 0,1120 = — 0,1120. 
er dy/g=25° x dop)g=200 R 


Multiplizieren wir ihn mit > Ay = 0,04365, so erhalten wir Jf = — 0,0049. Da f (25°) = 0 


ist, so wird f (20°) = 0,0049. 
Die Werte, die wir eben berechnet haben, sind in der Tabelle 2 in. den weiteren Spalten 
eingeklammert eingetragen. Es sind dies die Werte für 


2 fe a) 1 ( d N af af e N) 
h*3 ‚. — AlhR#—), U 70 —, 21 BE FA df. 
dp do mittel” YA : do Ei dp’ do do mittel j an f 


1b, 
Bei der Berechnung von — (7 En ea of hatten wir den Betrag = EB: > 7*3 f (20°) nicht be- 


Pl mittel 
rücksichtigen können. Nun nn wir gr näherungsweise berechnen und die ganze Rechnung 
noch einmal wiederholen. 
*\3 
Für = 20° wird Re h*2 f (20°) = 2) 120°) = n 1,192 0,0049 = 0,0015. 
x bi Lan 4 4 

Die verbesserten Werte sind in der Tabelle 2 ohne Klammern eingetragen. Wie sich zeigt, hat 

af 


diese Verbesserung auf die Werte Pen und f keinen Einfluß. 


1 
Nun gehen wir zur Berechnung der Werte für = 15°. Den Wert für (z nr f *) r 
. g= ° 


schätzen wir zu 0,0060, um die Berechnung gleich genauer durchzuführen. Die Berechnung gab 
den Wert 0,0053. In der Tabelle 2 sind weiter nur die endgültigen Werte eingetragen. Die Be- 
‚rechnung wird soweit fortgesetzt, bis wir zum Werte f(@,) = 0 gelangen. Für unser Beispiel er- 
‚halten wir auf diese Weise 
E- 2 = — 75°, 


* 
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8860 JS 


EE010 


62010 


6860°0 
1260°0 
LE80°0 


‚zrL0°0 
2590°0 
.7C0°0 
8500 
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D. Größe und Richtung der Lagerkratt 
Nachdem wir durch numerische Berechnung für unser Beispiel die Funktion f (P) ) und ihre 


Ableitung un gefunden haben, wollen wir Größe und Richtung der Lagerkraft bestimmen. 


Wir finden zuerst nach Gleichung (8) die Werte a, a, und a;. Zur Prüfung finden wir a, 
noch nach einem zweiten Verfahren. Es war 


> a h*f’do. 

Hieraus folgt durch ne. | 

E" | | o=— [hr fdp. : 

Weiter war 9, x | | 

a a 
P 


und 9 
d, — [ Beute do . 
Po 
Die Integranden sind im Bilde 3 dargestellt. 
} 
1d are 02 I 
ES N a SURT, SI TE 
k | Ser, b 
2 RT RE Er Se 0 % 
Bila 3 | A 2 we 
‚Durch a Die erhalten wir 
ee = 0,0508 nach dem ersten und _ | Ya 
2 ER TR = 0,0517 nach dem zweiten Verfahren, | nr 
Say im Mittel a, = 0,05125. } jr 
eg Weiter wwde 4, = 0,0435 und Bu 00080. 15. ‚ 
EN H ÄE b | | 
Br Re S Nunmehr ist Ban Gl. (18) mit Fe = 1: 
TER N "a 
! GEHT: 
. REP, ber R r ö an i % « 
K _% un — ist, so folgt rn u: Re 
. a. gr > da- ‚ 1 1 a 
Be r = a, +42 Be, 
u Br Beispiel ist diese Bedingung trotz der vielen Operationen EL gut erfüllt. Br 


igi Issen; gie Zahlenwerte @g, 4, und a, und wählen BR 3 
eR- EM e 00514, =, 0434 und a, = 0,0080. | 


Re: 
N] 


= = 0181 und = 2,33. 


190 Weber, 


Da «a zu 2,33 ea war, so wird B= 


Nun bestimmen wir Größe und Richtung = Lagerkraft. Da die Näherungsfunktion 


' | p* ein Produkt von f(p) und g (y) ist, ist die Bestimmung recht einfach. Wir benötigen zur Be- 
Br, rechnung das bestimmte Integral ' 
+B Br i 
so | 
2 op - P—XTgap). 
a 
Für unser Beispiel wird 
1,184 
nn 184 (3,77 _ Tg 3,77) = 2,83. 


2.3500 


Nun finden wir die Kraftkomponenten Pi in Richtung = 0 und P, in Do = - 900 
und erhalten mit Gl. = ) 


one Ba 
Po= p cos p d2z da = ——— Ip) g(y) cosp dp dy = 


ER 


AN) 
4 Ba 
WAT ERE bo Jg) cosp dp 
P% : 
und ebenso 
(3) Q 
6nu, 5 ’ 
Dr Ra 


ICYe 
Bild 3 zeigt die beiden Integranden für unser Beispiel. Die numerische Berechnung ereaD: 


Ti ) cospdp = 0,7830 und ro p) sing do = — 0,1035 . 


Mit 7(0,1830% + (0,1035% = 0,210 und gp=are ei = — 0,5656. — — 32,410 Dar 


erhalten wir die Lagerkraft 


a) 


So 3 Bi 
P= —— 2,83 0,210. 


Der Winkel zwischen der Kraftrichtung und der Richtung = 0 ist sgeseben durch den be 
rechneten Winkel 9p = — 32,41, 


,\ ; Ce 
P 


En: 


Führt man die Berechnung für eine genügend große Anzahl von Werten % Pr. "und Ha) ‚durch, r 
so erhält man ein vollständiges Bild für alle auftretenden Verhältnisse. % DER a E 


BE WIE NE für die Funktion f(p) bei vorgegebenen 5 und. 2 die Banden ge 

S(p)= 0 und — A Se - =( an. Hierfür en wir den dazugehörigen \ 
Sind die ee 9, und 1 gegeben, für die (op) = 9 und fi (9 

(91) 


BAnEHSON BER CZ 


—=0 fort. Dann wählen wir wieder wie ie vorher 
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Ein Verschiebungsplan für das räumliche Fach- 
werk. 


Mit Hilfeder dualen Abbildung werden die in räumlich 
gelegenen Ebenen auftretenden Verschiebungsvorgänge 
der Raumknoten auf ebene Probleme zurückgeführt. Da- 
durch wird es möglich, in zwei ebenen Bildern (Grundriß 
und Dualbild) den räumlichen Verschiebungsplan eines 
Fachwerks im Raum rein geometrisch durchzuführen. 
Das Verfahren hat zahlreiche Zeichenkontrollen, die der 
exakten Durchführung gute Stützen und dadurch eine 
ziemliche Sicherheit geben. 


Analog dem Williotschen Verschiebungsplan für das 
ebene Fachwerk läßt sich auch für das räumliche Fach- 
werk ein Verschiebungsplan aufstellen. Die zeichne- 
rischen Schwierigkeiten, die einer graphischen Dar- 
stellung in Grund- und Aufriß entgegenstehen, sind 
hauptsächlich darin begründet, daß die Ebenen, in 
denen sich die Verformung abspielt, nicht zeichen- 
gerecht liegen, weshalb sich die Behandlung durch Um- 
projektionen bzw. Umklappungen so umständlich ge- 
staltet. Die seitherigen Lösungsversuche, dieses Raum- 
problems bedienen sich zur Erreichung brauchbarer 
Methoden eines Abbildungsverfahrens des Raumes in 
die Ebene. Die zugeordnete Abbildung von B. Mayor 
und R. v. Mises!) und die von W. Prager?) waren 
Ausgangspunkte für zwei Lösungen der Grundaufgabe 
der kinematischen Knotenpunktsverschiebung des 
Baumfachwerks°®), die jedoch noch recht unüber- 
sichtlich sind, da die Ergebnisse im wesentlichen mit 
Doppelabbildungen von Geraden gewonnen werden. 
Mit den Grundlagen der dualen Abbildung von 
R. Sauer*) lassen sich nun die räumlichen Vorgänge 
rein geometrisch in die Abbildungsebene übertragen und 
eine klare und übersichtliche Konstruktionsgrundlage 
für das Verschiebungsbild gewinnen. Die Zugänglich- 
keit der räumlich orientierten Ebenen wird dadurch 
erreicht, daß sich jede Ebene als Punkt abbildet und die 


Schnittgerade zweier Ebenen als Verbindungsgerade . 


der zugehörigen Punkte erscheint. Mit Hilfe dieser 
Abbildungsgeraden, die den entsprechenden Geraden 
im Ausgangsbild eindeutig zugeordnet sind, werden die 
Verschiebungen im Raume in einem ebenen Bild durch- 
führbar und liefern Ergebnisse, die ohne Schwierigkeit 
wiederin das natürliche Bild übertragen werden können. 


1. Die Abbildungsgesetze °) 


Die duale Abbildung benutzt die Darstellung eines 
Vektors im Raum in eine Strecke einer festgelegten 
Abbildungsebene. Nennen wir die Abbildungsebene x, y 
(sie liege z. B. parallel zur Grundrißtafel des Raum- 
bildes), dann bilden sich ab 

Ser yayrıkeamcızt, 
x,y sind die Koordinaten bzw. Komponenten der 
Strecke lin der dualen Abbildung, /’ ist die Länge der 
Streckenprojektion auf die x,y Ebene; «*, y*, z* sind 
die Komponenten des zugehörigen Raumvektors, c ist 
eine den Zeichenverhältnissen anzupassende Abbildungs- 
konstante. a ist dann der Abstand des dualen Bildes 
von einem festgelegten Abbildungspunkt A. Geome- 
trisch gesehen werden also die z*-Komponenten räum- 
licher Vektoren in Abstände a von diesem Abbildungs- 
punkt verwandelt. 
ey 2* 


IT. m=Ue Tm—— 


Var pr 


Da es hierbei nicht auf die absolute Länge der Vektor- 


strecke ankommt, erscheint die duale Abbildung einer 
im Raum orientierten Geraden als Bild einer Geraden 


.*) R. v, Mises, Z, Math. u, Phys., Bd.64 (1916), 5.209. 

2) W, Prager, Z.angew, Math. Mech. Bd.6 (1926), 8.341. 
3) A, W. Sotoff, Z. angew, Math, Mech., Bd,8 (1928), 8,395 
und F, Reinitzhuber, Ing.-Arch, VI (1935), 8.318. 

*) R. Sauer, Z, angew,. Math, Mech., Bd.20 (1940), 5.174, 
» ) Vgl, auch H. Dietz, Ing.-Arch, XVI (1947), 8.13, 


in der Abbildungsebene. Beide Bilder sind eindeutig 
einander zugeordnet. Eine Ebene, die sich im Raum 
durch zwei Geraden kennzeichnen läßt, hat in der 
Abbildung dann den Schnittpunkt der beiden zuge- 
ordneten Abbildungsgeraden als zugeordnetes Bild, und 
es entsteht ein ganz allgemeines Abbildungsgesetz, in 
dem 
die Geraden des Raumbildes 
zu Geraden des Dualbildes, 
die Ebenen des Raumbildes . 
zu Punkten des Dualbildes 
und die Punkte des Raumbildes 
zu Ebenen des Dualbildes 
werden. Die Ebene erscheint in der Abbildung als 
Punkt mit den Koordinaten b;, b, (d.h. im Abstand 


b= Vs; -+ by vom Punkt A), die definiert sind: 


ey ex Fi a at 
en A ee a A NeE 2 
by = —e e, ? by Fe 6,’ b e; Ve + ey: 


Darin ist die Ebene (oder eine begrenzte Fläche in 
derselben) mit e und ihre Projektionen in der x-, y- 
bzw. z-Richtung sind mit e,, ey, e; ‚bezeichnet. 

Eine Ebene e senkrecht zu einer Geraden g bildet 
sich als Punkt ab, der auf dem Lot zur Geraden durch 
Punkt A liegt und dessen Entfernung b vom Abbil- 
dungspunkt A zu der Entfernung a der Geraden g 
von A im geometrischen Verhältnis steht: 


Bild1. Duale Abbildung einer Ebene e, 
die senkrecht zur Geraden g steht 


Benennen wir in Bild 1 die Gerade g und die Kompo- 
nenten einer beliebigen Vektorlänge auf g mit gx, 9y, 9z- 
Die Ebene e mit ihren Projektionen e;, ey, e; stehe 
senkrecht auf der Geraden g. Dann ist nach den 
Abbildungsgesetzen 


g2 . rn 
ac: „it ’=Vg2+g 
%y Ww%& I 


er also, wird 
[23 g2  % 92 


und wegen der Richtung 


; a 1 “: Sl 9x, 
u rt ER 
; 
n 9 
Fa Er : 
b= ra Ve; + DE und damit 
et BE 0 
a.b= di bzw. 7 =y° 


Diese Erkenntnis hilft zu einer vereinfachten Kon- 
struktion des Abbildungspunktes für eine zu einer 


‘Geraden senkrechten Ebene und, da die Dualbilder für 


Gerade und Ebene rein graphisch gefunden werden 
können®), läßt sich auch die ganze Konstruktion des 
Dualbildes und des folgenden Verschiebungsplanes auf 
rein graphischem Weg durchführen. 


*) vgl. W,Schlink, Techn, Statik, 4. u,5, Aufl, Springer 1948. 


a Kal 


, 
= 


E 


; 


u > 
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2, Das Verschiebungsbild des Dreibocks 
(einfacher Raumknoten). 


Die Grundlagen des Verschiebungsplans seien zu- 
nächst am einfachen Dreibock erklärt. Jeder Stab hat 
infolge seiner Stabkraft eine Längenänderung 

AEWSHTL 
A=yF 

Die Forderung, daß der Knotenpunkt auch nach der 
Deformation wieder alle Stäbe vereinigt, verlangt, daß 
jeder verkürzte oder verlängerte Stab sich mit seinem 
oberen Ende auf einer Kugelfläche — mit dem Stab- 
fußpunkt als Mittelpunkt — bewegt. Im Schnittpunkt 
der drei entstehenden Kugelflächen finden sich dann 
die drei Stabenden wieder zusammen. Die Kugel- 
flächen lassen sich wegen der Kleinheit der Ausschläge 
und der verhältnismäßig großen Radien als Ebenen 
senkrecht zum Stab darstellen. Die Ebenen treffen sich 
im neuen Knotenpunkt 0*. Operieren wir nun mit 
diesen Ebenen und ihren Schnittgeraden mit den 
Ebenen der Stäbe, die sich im Dualbild alle als Punkte 
bzw. als Verbindungslinien dieser Punkte darstellen, so 
kommen wir auf bequemem Weg in einem ebenen 
Lösungsbild zur Verschiebung s des Knotenpunktes. 
Ebenso finden wir leicht den senkrecht zur Stabachse 
stehenden Weg des einzelnen Stabendes. 

In Bild 2, das die Spitze eines Dreibocks zeigt, werden 
zunächst die Fußpunkte der Stäbe ohne Verschiebung 
angenommen. Damit wird die Lage des neuen Knoten- 
punktes 0* nur von den Längenänderungen der Stäbe 
ahhengig, die wiederum durch ihre Stabkräfte bedingt 
sind. 

Die Ebene I, die senkrecht zu Stab ® steht, schneidet 
z. B. die Ebenen e (aus @und ®) und f(aus @und ®) 
in je einer Geraden 1e und 1f. Ebenso lassen sich die 


Ebenen II und III, senkrecht zu den Stäben @und ®& 
darstellen durch ihre Schnittgeraden 2e, 29 bzw. 3/ 
und 39 mit den angrenzenden Stabebenen f und g. 
Gehen wir nun von den Schnittpunkten dieser Geraden 
senkrecht zu den zugehörigen Stabebenen, so müssen 
sich diese Lote im neuen Knotenpunkt 0* des ver- 
formten Dreibocks schneiden. Kinematisch gesehen 


lassen wir bei diesem Bestimmungsvorgang erst sich 


die Ebenen je zweier Stäbe verformen und schwenken 


_ diese verformten Ebenen aus nach dem Treffpunkt 0*, 
Daß diese Operation mit der Bewegung des verschh- 


benen oberen Stabendes auf einer Ebene senkrecht zum 
Stab (als angenäherte Kugelfläche) identisch ist, läßt 


sich geometrisch leicht überblicken, da alle Vorgänge 


in diesen senkrechten Ebenen liegen. 

Im Dualbild, das uns die Durchführung dieses räum- 
lichen Verformungsbildes in einer ebenen Darstellung 
erlaubt, zeichnen wir uns zunächst die duale Abbildung 
der Stäbe und der darauf senkrecht stehenden Ebenen f 
nach den oben erwähnten Abbildungsgesetzen. Als | 
Abbildungsebene wählen wir aus Einfachheitsgründen 
den Grundriß (die Lage der Abbildungsebene ist be- 
langlos für den Lösungsgang). Die Geraden le und If 
sind die Schnittgeraden der Ebenen e mit I bzw. f 
mit I, erscheinen also jetzt im Dualbild als Verbindungs- 
geraden der entsprechenden Ekenenpunkte. Die Ver- 
schiebungen der Stäbe in den Stabebenen können nun 
mit den Parallelen zu diesen Geraden eingezeichnet 
werden und liefern die Schnittpunkte Z, F und G in 
den Ebenen. Das Lot im Punkt Z auf der Ebene e . 
ist zugleich in der Ebene senkrecht zu Stab @undin j 
der Ebene senkrecht zu Stab ® enthalten, entspricht 
somit der Verbindungsgeraden I—Ilim Dualbild. Diese 
Lote für die drei Ebenen treffen sich im Grundriß dann 
im neuen Knotenpunkt 0*. (Die dritte Gerade stellt 
schon eine Zeichenkontrolle dar!) 

Aus Dualbild und Grundrißfigur lassen sich ein- 


- deutigdiezugeordnetenräumlichen Strecken bestimmen. 


Wir erhalten allgemein für die z*-Komponente einer 
bestimmten Strecke I: 
2% ’ 
N 
6 
worin a dem Dualbild und /” dem Grundrißbild des 
Systems zu entneh- 
men ist. (Die Aufriß- 
figur wird dabei 
' überflüssig, sie ist 
für den Lösungsgang 
durch das Dualbild 
ersetzt.) Die wirk- 
liche Größe der 
Strecke wird 


tr Yı+(e). 


Damit lassen sich 
nun die für statische 
Untersuchungen in- 
teressierenden o- 
Werte, d.h. die Be- 
wegungen des Stab- 
endes senkrecht zur 
Stabachse (in der 
Schwenkebene), und 
der s-Wert, d.i. die 
Verschiebung des 
Knotens, in ihrer 
wahren Größe be- 
stimmen aus den 
Daten der beiden’Bil- 
der (Grundriß un 
Dualbldl) 


Be RN OR. 
i T e ' Yı-@) A 6 


m. 


Stäben enthalten sind), auf den drei 
Punkte anschneiden die 1. auf Er 
müssen und deren Verbind yaralle 
verschiebung s’ im Grundri 

schnittenen Punkte sin d die, 
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zelnen Stäbe, die sich alle drei in der gemeinsamen 
Geraden s treffen. 


.* 3. Der Knotenpunkt im Raumfachwerk 


Beim Knotenpunkt in einem Raumfachwerk sind 
die Fußpunkte der Stäbe an anderen Knotenpunkten 
angeschlossen, die Stäbe erleiden also schon von diesen 
her Verschiebungen. Wenn wir von diesen Fußpunkten 
im Verschiebungsplan aus nun die neuen Stabdeh- 
nungen antragen, gibt das Verschiebungsbild des folgen- 
den Knotens keine geschlossene Figur, d.h. die Lote 
auf die Stabebenen schneiden sich nicht mehr in einem 
Punkt. Das Antragen der Stablängenänderung an die 
Fußpunktslagen mit den unveränderten Stabrichtungen 
stellt eine Parallelverschiebung der Stäbe dar, die, von 
den neuen Fußpunktslagen ausgehend, sich jetzt nicht 
mehr im oberen Knotenpunkt schneiden, sondern wind- 


g 
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als die statisch interessierende Größe für den Stabaus- 
lenkwinkel von Bedeutung ist. o steht senkrecht auf 
der Stabrichtung und liegt in der Spur der Ebenen aus 
Verschiebung des Fußpunktes und Stab einerseits und 
der zum Stab senkrechten Ebene andrerseits. Dement- 
sprechend finden wir die zugehörige Richtung als Ver- 
bindungsstrecke der beiden Ebenenpunkteim Dualbild, 

Mit der Lösung des Verschiebungsecks erhalten wir 
für jeden einzelnen, Stab sowohl die o-Werte am Fuß- 
ende als auch an der Spitze, die für die Feststellung der 


Stabdrehwinkel y = E von Bedeutung sind. Weiterhin 


liefert das Dualbild die Lage der Verformungsebene als 
Schnittpunkt der Knotenverschiebung s mit den Dual- 
bildern der Stäbe. Die wirkliche Lage im Raum läßt 
sich für oe sowohl wie für die Ebenen leicht nach den 
Gesetzmäßigkeiten der Abbildung hestimmen. 


Bild. 3 Fachwerkbeispiel in Grund- und Aufriß 


schief aneinander vorbeilaufen. Um nun 
hier einen Lösungsweg zu erhalten, ohne 
daß wir neue Stabrichtungen bestimmen 
müssen (die Richtungsänderungen | 
der Stäbe sind von höherer Ordnung 
klein und erscheinen nur in dem 


stark vergrößerten Verschiebungs- 0: __w __w __ 


57 Re 
m 


mn 


c«5cm 


eck von Bedeutung), legen wir die 
. Stabdehnungsrichtungen der Stäbe 
stets vom Ausgangspunkt 0 aus- 
ehend, nehmen den senkrecht zur 


i Kräftemaßstab 


 Stabrichtung liegenden Anteil der 
Fußpunktsverschiebung weg und tra- 
gen von diesem so gewonnenen Punkt 
aus die Stabdehnungen an. Das Ver- 
 schiebungsbild des neuen Knotens 
wird dann in der beschriebenen 


Weise durchgeführt. Diese Maßnah- 


menentsprechen imräumlichen Fach- 
werk einem Zusammenfassen der drei 
Stäbe im oberen Knotenpunkt mit 
_ einer Verschiebung der Fußpunkte um 
‚das Schwenkmaß o für den Stab, das 


te, 


Bild 4. Graphische Bestimmung der Stabkräfte mit der dualen Abbildung 


* 


Abbildungskonstante 


4. Beispiel 


(Bild 3), das außer seinem Bildungsgesetz keinerlei 
Sonderbedingungen unterworfen ist, ist die Stabkraft- 


bestimmung mit dem Dualbild durchgeführt (Bild 4). 


Der einzelne Raumknoten mit den drei unbekannten 
Stabkräften durch einen Punkt wird dabei zum ebenen 
Problem dreier in der Abbildungsebene zerstreutliegen- 
der Bildkräfte. Die Gleichgewichtsbedingungen über- 
tragen sich entsprechend ins Dualbild?) und erlauben 


2) Siehe H. Dietz a.a, 0, 


E 


Bild5. Verschiebungsplan der F 
« ; 


a 


y 


Für ein einfaches Raumfachwerk mit 6 Knoten 


folgende Tabelle: . N 


RE 
chwerkkart 


eine graphische Lösung der Stabkräfte, deren Ergebnis 
die Grundrißprojektionen / der Stabkräfte sind. Für — 
die Längenänderungen der Stäbe im Grundrißbild 
brauchen wir nicht auf die wahren Größen der Stab- 


kräfte zurückzugehen, denn es ist 
’ 


(02) 


.] % So 1 ‘ } R I“ ; | 
a BF’ also kA = 7. 7 (oder 7 a 
bei errechneten Stabkräften). Bi Bi, 
Die Abmessungen und Werte für die Stäbe, die dem 


Verschiebungsplan zugrunde gelegt sind, enthält die 


! 
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Absolute Längen in cm der 


Zeichnung (Kräfte kg) Ergebnisse 


Stab 
(Kraft) 


PAST: 


10°3cm 
11583. 07802178 8 — 5,580 
2| 2| 3 | 4 | 3,606| 5,385 — 0,946 
31 2| 3|4| 3,606| 5,385 — 0,946 
4110| 3 | 4 110,44 |11,18 —+ 4,880 
5I 2/3 | 4 | 3,606] 5,385 — 1,770 
6] 2| 3 | 4 | 3,606) 5,385 — 0,942 
Ze0 6 1 0:1.6 6 — 0,390 
8112| 0|0|12 12 — 0,540 
9] 2| 3 | 7 | 3,606] 7,874 — 2,008 
10/1121 0|0|!12 12 — 9,756 
18120 6 | 04,76 6 —124 | — 0,744 
12] 2| 3 | 7.| 3,606] 7,874 —170 | — 1,340 
13] 2) 9 7 | 9,22 111,576 + 47| + 0,544 
14| 2| 3 | 7 | 3,606) 7,874 —181| — 1,425 
15] 14| 3 | 7 |14,32 |15,94 —357 | — 5,690 
16| 2| 9 | 7:|° 9,22 |11,576 + 43| + 0,498 
17) 2| 3 | 7| 3,606] 7,874 —541 | — 4,260 
18] 14| 3 | 7 |14,32 [15,94 —+720| +11,470 
P 500) 0 1400! 500 640,3 14 
P,;,| 0 0 4001 0 400 M =400 5 cmkg 
P,| 0 0 3000 0 1300 M = 300 :5 cemkg 
P,| 0| 0 300 0 300 M = 300 - 5 cemkg 
E 3,75 


3,75 

Dabei ist für alle Stäbe der gleiche Querschnitt an- 
genommen und das Ergebnis 8’ -I in einem zeichen- 
gerechten Maßstab verzerrt, der nach der Dimensio- 
nierung leicht auf wirklichkeitsgerechte Werte umge- 
rechnet werden kann. ei 

Die Verschiebungen der Fachwerksknotenpunkte 
a—f sind nun, der Reihenfolge des Fachwerkaufbaus 
entsprechend, nach den oben gegebenen Regeln kon- 
struiert und ergeben den in Bild 5 gezeigten Verschie- 
bungsplan. Die Knotenverschiebungen s sind im 
Grundrißbild als Verbindungsstrecken s’ der Knoten- 
punktslagen a—f mit dem Nullpunkt 0 abzugreifen, 
die 0 -Werte findet man als Strecken zwischen den 
Endwerten von A’ einerseits und den verschobenen 
Knotenpunkten a—/f andrerseits (für Knoten b sind die 
Werte eingezeichnet). Ausder Streckenlänge im Grund- 
riß und den Abständen von A im zugehörigen Dualbild 
kann man sichleichtdie wahre Größe der Verschiebungs- 
werte und ihre Richtungen, bzw. ihre Ebenenerrechnen. 

Darmstadt. H. Dietz. 


Bemerkung zur Fourierschen Reihe für die 
dreidimensionale Greensche Funktion der 
Wellengleichung des unbegrenzten Raums 
in den Koordinaten des Drehparabols. 

In einer älteren Arbeit des Verfassers!) wurde die 
mathematische Beziehung für die Greensche Funktion 
exp(ikR)/(ikR) der Wellengleichung in den Koordi- 
naten des Drehparabols bei beliebiger Lage von Quell- 
punkt und Aufpunkt mit dem gegenseitigen Abstand R 
in der Form der Gl.(1, 1a, b) hergeleitet. Darin 
bedeuten %k die reelle oder komplexe Wellenzahl, 


ER, NT cosp (PH) 
ikR rer. pip! 
—0y io 
; ER 2) . (1) 
ri n(s+ 2 
{ —0y— ix 
p+1 Kerr 
rn (- 04252) Panee, nf, m) de 
l+p» 
(al) 


—I7H. Buchholz, Zeitschr. f. Phys., 124 (1948), 8.196 
bin218. 
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Ps, 212 (&, 7, &o, N) RAN 
& >$) . 
uf 2 . fs 
M+s,212(—iE)  W+snl2(— 8 &) 
ie Ti 


| 
Mi» (—i&) Wisoe(—i®) 
LE 
(1a 


(1 8) , 
(& <$) 
Mm>m , 
M_-s.»12 (— in) W-s.n2 (— im) 
im im | 
"IM-sm2(— in) W-sr2 (— in‘) | 
in)"; if): 


(N <n) 


) 


&, n, g oder &,, 70; 9. die parabolischen Koordinaten 
von Aufpunkt und Quellpunkt, die mit den recht- 
winkligen Koordinaten x, y, z oder den Zylinder- 
koordinaten 0, p, z oder den Kugelkoordinaten r, ®, p 
über die Gl. (2a, b,c) zusammenhängen, und Ms,»/2 (2), 
Ws,»/2(2) die beiden konfluenten hypergeometrischen 
Funktionen von E. T. Whittaker, die wir hier auch 
als parabolische Funktionen bezeichen werden. Gemäß 
den Angaben zu Gl.(la) setzen sich z.B. im Falle 


2=0-'cosp=rsin®#-cosp—2YEn-cosp - (2a) 
. (2b) 
(2e) 


eines &>E, 19» <n die Kernfunktionen Ps,p/2 im 
Integranden von (1) jeweils aus dem Produkt der vier 
parabolischen Funktionen Ms,»/2 (-i&), Ws,»/2(-i &5) 
M_s,p/2 (-in4) und W_s,p/2 (-in) zusammen. Ist 
&,<E, und 7,< n, so hat man die Funktion Ps,p/2 
aus dem unteren Produkt der ersten und der zweiten 
geschweiften Klammer zu bilden. Damit dürfte die 
Schreibweise der Gl. (1a) in ausreichendem Umfange 
verständlich gemacht worden sein. Eine Umformung 
der Integrale von (1) in unendliche Reihen ist bei 
einer Zusammensetzung des Integranden, wie sie z. B. 
dem Falle &,>&, 7 >n entspricht, nur möglich, 
wenn entweder 


VEE>VE+Ym+Yn oder Ym>VEo+VE+ In 


ist, je nachdem nämlich der Integrationsweg nach 
rechts an die positiv reelle Halbachse oder nach links 
an die negativ reelle Halbachse herangebogen wird. 
Wir behalten daher die in allen Größenfällen von 
&,n...gültige Integraldarstellung bei. 

Für die Anwendungen auf Physik und Technik ist 
esin der Regel vorzuziehen, was ja auch A. Sommer- 
feld stets befürwortet, die in (1) auftretende Fourier- 
Beihe nicht mit der Funktion cos anzuschreiben, son- 
dern durch Zählung der Summationsvariablen p von 
— 00 »*+-+00 mit der Exponentialfunktion, weil da- 
durch die lästige Ausnahmestellung des Gliedes p = 0 
aufgehoben wird und obendrein eventuelle Rechen- 
operationen an der Variablen p leichter ausgeführt 
werden können. Es ist in der Hauptsache das Ziel der 
vorliegenden Arbeit, die Schreibweise von (1) dieser 
Forderung anzupassen. 

Denkt man sich im Sinne dieser Aufgabe die Funktion 
cos (p [P—9u]) in die halbe Summe der beiden Exponen- 
tialfunktionen exp(+iplp— 9) und exp(—ip 
[p — 9]) zerlegt, so hat man es darnach auf der rechten 
Seite der Gl. (1) nach dem durchweg vorzunehmenden 
Ersatz von 1 + öop durch 2 zunächst mit einem ein- 
zelnen Integral zu tun, das dem Wert p = O entspricht, 
und außerdem mit zwei unendlichen Reihen, in denen 


y=o-sinp=rsind-sinp=2Y&n-sinp 


ss =2=r-. cd =E—n 


Fangen 


Kr} 


Be u ine FU un Be» en 
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jedesmal p von 0... zählt, von denen aber die eine 
vor dem Int egralzeichen den Faktor exp(+ip[p— pol) 
enthält, die “andere den F aktor exp(—ip[lP—@o]). 

Um diese drei Bestandteile zu einer einzigen Reihe 
zusammenzuziehen können, in der p von — über 0 
bis oo zählt, müßte es möglich sein, in der Reihe mit 
dem Faktor exp (—ip- [pP —9,]) mit p=1-: +0 


von + p zu — p überzugehen, so daß nunmehr p von 


—1-::—o zählt. 

Von diesem Übergang von + p zu — p werden dann 
aber auch die Faktoren in den einzelnen Integranden 
betroffen. Bei der Funktion W,./2 (z) hat das nichts 
auf sich, denn sie ist für alle Werte u eine in « gerade 


Funktion. Bei der Funktion Mz,/2 (2) ist aber dieser 


Übergang gemäß der Definitionsgl. (3) dieser Funktion 


Ii+u 
Mr,u2(@)=z ? -e"?.,P, Br Li :) 
1+u & (+) BR 
_y* Park EN 8) 
Be; (l+Fu),:A! 
@=1 (y=ala+l)—(a+4—1) 
"nicht ohne weiteres möglich, denn für ein u= —p 


würde ihr zufolge im Nenner der Kummerschen 
Funktion ,F, eine negative ganze Zahl zu stehen 
kommen ade daher würden in der unendlichen Reihe 
alle Glieder mit A > p unendlich werden. Nun ist aber 
wegen 1/I(—p)=0 für p=1,2,3... 


De (8; 1+ß; 2) 


— Jim sa Ze FA ale 
N — () 


(a) T(L+B+R) 2 


B I‘& + )) zi 
— Jim ] E i 
30 Pe Po) PT(l—p—d+N 4 


e ) 
wa ToTa—p—5+n Mf 


Bau NT DIET PA EN BC) 
Ton tpHl) A Te) 


zP 
uhetmp+ln. 


Mithin ist im Hinblick auf die Gl. (3): 


149 ) 
Tr 
lim Ms, u/2 (z) .y ( 2 7 \ Mz, »12(2) 5 


„-— „I (l+M) De ) P(1+p) 


- (4) 


(p=1,2,3.--») 


. Bei dieser Sachlage empfiehlt es sich, an Stelle der 
in (1) und (la) auftretenden Funktionen Mx ‚u/2 (2) 
und W,„/2(z) in dem hier gegebenen Zusammenhang 
mit den durch die Gl. (5a, b) definierten Funktionen 

m“ (z) und w) (2) zu rechnen. In Rücksicht auf die 


eben besprochenen Eigenschaften der beiden Whitt- 
akerschen Funktionen gelten dann für sie die Formeln 


Win — le . Mui2l2) 
m." (2) =2 Ta (5a) 
Wu 2... (bb 
r(Zr— )m Na x) m\* (2) 
?=01,2:»») (56) 


wir) Mut (yibel.) ER AR vl, 
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(5%, ß). Der Übergang } von + p zu — p ist Be ine ; 
dem hend umgeänderten Integral von (1) ohne 


1) 


weiteres durchführbar, und es entsteht damit für die 
Greensche Funktion 'exp(ikR)[(ikR) die allgemeine 
Darstellung (6) mit der Kernfunktion @5,7/2 von 1. Gl. (6a). 


ir (e-M) 


gikR SE 
fikR= —2: y: u 


Pp=-x 


_0,+i» 
[ehr | 


.- io 
27 [ 


Qs,pl2lE, 7; 8,10): ds 
1 = an 
(o “in ARE sn) 


Qs,n2 (€, Ms ;£, 70) 


(& (6a) 
(7 < m) a 
Bent (m) (in, 
E jm- 2 in Me w i 70) \ 
Ki m im)-w en] 
>) 


Ergänzt En 5 G1.(6) durch die ne Dar- 


stellung von e'"*jikR in Zylinderkoordianten, multi- 
pliziert die entstehende dreifache Gleichung mit 
exp(—ing) mit n = 0,1,2... und integriert sie über 
p zwischen den Grenzen 0. 2a, so gelangt man zu 
der Formel (7). In dem darin vorkommenden Integral 


- mit den beiden Zylinderfunktionen weicht der Inte- 


grationsweg auf der Seite 6<0 dem Punktt= —1 


2r \ + 


ikR : 
e er dpy—n: „ginn erklemzol-t 
ikR 

0 n x RT 


(@> 0) 
E (ikon: EZ). . m (iko ya=ı) N. 
. di. 
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wenn man darin 0,—>0 gehen läßt. Wegen 9 = 2 
- VEu no entspricht dem entweder der ‘Grenzübergang 


&,—>0 oder 9 >0. Für &,— 0 ergibt sich z. B. 


n -+2 BEZ? 


EURO er. 
ee FE 4 
(mr "9 
601,4 io 
1 ER 1-+n R 1l-+n 
2zi j n(s+ 2 en + 2 ) 
—-01,-io 3 


wi i8) m) (— in) w) (— im): ds (8). 
as . . KEANE 
nf et Pr re YR— 1) 


—nn 


(12. (n=0, +1, +2---) 


l+ \ 
(m> N5 lejnıl > Ze) 


n=0, +1 +2...) 


Eine weitere Transformationsgleichung entsteht, 
wenn man für exp(ikR)/(ikR) die auf ein Kugel- 
koordinatensystem bezogene Integraldarstellung heran- 
A Hier ergibt sich dann auf dieselbe Weise wie 
oben 
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die G1. (9), zu der man noch die letzte Zeile von Gl. (7) 
hinzuzufügen hat, um die eigentliche Transformations- 
gleichung zu bekommen. In (9) ist €’ ein Integrations- 
weg, der aus dem vierten Quadranten kommt, die reelle 
t-Achse zwischen den Punkten —1- - - 0 schneidet und 
innerhalb des ersten Quadranten wieder ins Unendliche 
ausläuft. Die P}”" sind die zugeordneten Kugelfunk- 


tionen. 


Darmstadt. H. Buchholz 


BUCHBESPRECHUNGEN 


Dr. Heinrich Blasius, Mechanik, Physika- 
lische Grundlagen vom technischen 
Standpunkt; zweiter Teil: Elastizi- 
tät und Festigkeit. VIllu. 230 S. m. 194 Abh. 
u. 64 Aufgaben. 3. Aufl.; unveränderter Abdruck der 
zweiten neubearbeiteten Auflage. Hamburg 1949. 
Verlag Boysen & Maasch. Preis brosch. 9,50 DM, 

In der vorliegenden dritten Auflage sind gegenüber 
der zweiten Auflage, welche in dieser Zeitschrift, Bd. 14 
(1934), S. 124/25 besprochen wurde, nur einige Druck- 
fehler verbessert. » 


Dresden. Neuber. 


Dr.-Ing. Max Schuler (beam. ao. Prof. a. d. Univ. 
Göttingen), Mechanische‘ Schwingungs- 
lehre. Teil I: Einfache Schwinger (Bücher der Mathe- 
matik und Naturwissenschaften, herausgeg. von Dr. 
Henry Poltz). 168 S. m. 123 Abb. Wolfenbüttel und 
Hannover 1949. Wolfenbüttler Verlagsanstalt. Preis 
kart. 6,— DM.-West. 

Auf 162 Seiten (in Format A 5) behandelt der Ver- 
fasser in drei Kapiteln die freien (ungedämpften und 
gedämpften) sowie die erzwungenen Schwingungen der 
Systeme von einem Freiheitsgrad. 

Im Kapitel über die freien ungedämpften Schwin- 
gungen wird ein besonderes Augenmerk dem Pendel so- 
wie anderen Vorrichtungen zur Zeitmessung zugewandt. 
Die Paragraphen über das physikalische Pendel und 
über die elastischen Drehschwingungen enthalten nicht 
nur Abschnitte über den Antrieb der Pendel- und Un- 


 ruh-Uhren, sie befassen sich auch mit der Quarzuhr 


und Vorschlägen zu einer „Lichtuhr‘‘. Dem (Schu- 
lerschen) Ausgleichspendel, dem Reversionspendel, 
dem „beschleunigungsfesten‘‘ Pendel ist u. a. jeweils ein 
‚Abschnitt, dem Zykloidenpendel und dem auf einer 


Schneide abrollenden Pendel jeweils ein Paragraph ge- 


zur 
Ei 


widmet, 

Im Kapitel über die freien gedämpften Schwingungen 
außer den üblichen Diskussionen des Ablaufs der 
n bei Anwesenheit von geschwindigkeits- 
onalen Dämpfungskräften, von Festreibung 


5 


und von quadratischer Dämpfung ausführlicher als 
üblich die Berechnung der Gleichgewichtslage einer 
Schwingung behandelt. 

Das Kapitel über die erzwungenen Schwingungen ist 
vergleichsweise überaus knapp und enthält nichts, was 
über das ganz Landläufige hinausginge. 

Die gesamte Darstellung ist: recht „schlicht“ und 
verwendet (abgesehen von einigen Besonderheiten) 
nur bescheidene Mittel; gelegentlich mangeltsogar etwas 
die Sorgfalt (im Rechnerischen und im Sprachlichen). 
Vor allem aber macht die Darstellung in keiner Weise 
Anspruch auf Systematik, sie verzichtet auch völlig 
auf die Aufdeckung etwas tieferliegender oder weiter 
ausgebreiteter Zusammenhänge. Diese Feststellung 
stimmt den Betrachter um so verwunderter, als ihm im 
Vorwort mitgeteilt wird, daß der Inhalt des Bänd- 
chens den ersten Teil des Stoffes einer Vorlesung bildet, 
die viele Jahre lang an der deutschen Universität mit 
der berühmtesten mathematisch-naturwissenschaft- 
lichen Fakultät abgehalten worden ist. 


Karlsruhe. K. Klotter. 


K, W,. Wagner (vorm. o. Prof. a. d. Techn. Hoch- 
schule Berlin), Einführung in die Lehre 
von den Schwingungen und Wellen. 
XVI u. 6408. mit 288.Abb. Wiesbaden 1947. 
Dieterichsche Verlagsbuchhandlung. Preis brosch. 
34,— DM. 

Von ganz anderer Qualität ist das zweite der hier 


‘zu besprechenden Werke über Schwingungslehre. Das 


Wagnersche Buch ist ebenfalls aus Vorlesungen 
entstanden; sie sind an der Technischen Hochschule 
Berlin abgehalten worden. Das Buch hat wie die 
Vorlesungen das Ziel, eine Einführung in das Gebiet 
zu geben. Was entstand, genügt aber erheblich weiter- 
gehenden Ansprüchen. Die Darlegungen bewegen 
sichauf einer Höhe, von der aus auch Zusammen- 
hänge aufgezeigt und Verbindungen zu Nachbargebie- 


‚ten geschlagen werden können. Das mathematische 


Rüstzeug, das gefordert: und geboten wird, ist diesen 
Absichten angemessen. 


Das 640 Seiten zählende Werk umfaßt einen bemer- 
kenswert reichen Stoff. Das erste Kapitel handelt von 
den Schwingungsformen, d.h. von der Kinematik der 
Schwingungen, die man ja zweckmäßigerweise an die 
Spitze stellt. Es enthält neben der Zusammensetzung 
und der Modulation von Schwingungen auch die Fou- 
rierschen Reihen und Integrale in dem weiterhin 
benötigten Umfang. Im zweiten Kapitel „Dynamik der 
Schwingungen‘, werden die freien und erzwungenen 


Schwingungen der Systeme von einem, von zwei oder 


von mehreren Freiheitsgraden untersucht; und zwar 
werden sowohl mechanische wie elektrische Gebilde als 
Beispiele herangezogen. Das dritte Kapitel handelt 
von den Schwingungen und Wellen in kontinuierlichen 
Systemen; von mechanischen Gebilden werden be- 
trachtet: Saiten, Stäbe, Balken, Luftsäulen, Mem- 
branen und Kugelstrahler, von elektrischen die homo- 
gene Leitung mit Eigenschwingungen und Wander- 
wellen, Potenzleitungen und Exponentialleitungen. Das 
vierte Kapitel ist einem speziellen Thema, den elektri- 
schen Wellen in der Ionosphäre gewidmet. Das fünfte 
Kapitel behandelt Schwingungsgebilde mit periodisch 
veränderlichen Elementen, das sechste Kapitel (recht 
ausführlich) die nichtlinearen Schwinger, ihre Eigen- 
schwingungen, ihre erzwungenen Schwingungen mit 
den charakteristischen Frequenzgängen und Instabili- 
täten und die damit erzielbare Frequenzvervielfachung, 
die Selbsterregung von Schwingungen durch Lichtbogen 
und durch rückgekoppelte Verstärkerröhren, die Kipp- 
schwingungen und die Erscheinung der Mitnahme und 
des Ziehens. Acht (ungleich lange) Anhänge bieten Er- 
gänzungen über Fouriersche Reihen und Integrale, 
die Grundregeln der Matrizenrechnung, die wichtigsten 
Regeln für das Rechnen mit der Laplace-Transtor- 
mation, über Klektronen-Plasma-Schwingungen, die 
Spiegelung und Brechung von elektrischen Wellen, so- 
wie über zwei kleine Gegenstände. 

Die Aufzählung mag einen Eindruck geben von der 
erwähnten Vielseitigkeit des Gebotenen, und es sei be- 
tont, daß die Probleme mit den jeweils angepaßten 
Hilfsmitteln angegriffen und bewältigt werden. Die 
Lehrerfahrungen des Verfassers zeigen sich vielerorts 
in vorteilhafter Weise. 

Wo viel Licht ist, werden aber auch, Schatten ge- 
worfen, Diese rühren z. T. von Eigenwilligkeiten des 
Verfassers her. So hält er sich z. B. nicht an die an- 
erkannten Definitionen über die Abzählung der Frei- 
heitsgrade, sondern schreibt z. B. einem Feder-Masse- 
System zwei Freiheitsgrade zu, weil seine Differential- 
gleichung von zweiter Ordnung ist. Auch die Benen- 
nung der Kopplungsarten (Feder-Kopplung, Massen- 
kopplung, kapazitive und induktive Kopplung) und 
ihre Behandlung lassen mit der wünschens- 
werten Deutlichkeit die Bedingtheit durch die jeweili- 
gen Koordinatensysteme erkennen. Es mag ferner be- 
fremden, daß im fünften Kapitel über die Gebilde mit 
periodisch veränderlichen Elementen die grundlegende 
Hillsche und Mathieusche Differentialgleichung mit 
keinem Wort erwähnt wird, so daß auch deren weit- 
entwickelte Lösungen und Lösungsmethoden nicht 
herangezogen werden können; die Differentialglei- 
chungen werden vielmehr sogleich mit der Methode der 
sukzessiven Approximation angegriffen, ’ 

Verwunderung muß dann noch erregen, daß gerade 
die dem Verfasser doch so überaus wichtige Frage der 
mechanischen und elektrischen Analogien in recht un- 
zulänglicher Weise behandelt ist. Es wird nur die land- 
läufige Zuordnung erwähnt, bei der die Masse der In- 
duktivität entspricht, nicht aber die andere, ebenfalls 
mögliche und überdies sehr viel zweckmäßigere, wo die 
Masse der Kapazität zugeordnet wird, obgleich Fire- 
stone schon 1932 deren Vorteile aufgezeigt hat.!) Zweck- 
mäßiger ist diese zweite Zuordnung, weil mit ihrer Be- 
nutzung bei Gebilden mit mehr als einem Freiheitsgrad 
eine Reihenschaltung wieder in eine Reihenschaltung, 
eine Parallelschaltung wieder in eine Parallelschaltung 


») F.A. Firestone, Journ, Acoust, Soc. (1933) 8. 249. 
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übergeht. Solange man bei Systemen von einem Frei- 
heitsgrad stehen bleibt, fällt diese Unterscheidung weg, 
dann ist das Problem aber nahezu trivial. Daß die Un- 
zweckmäßigkeit der benutzten (ersten) Zuordnung im, 
Buch nicht auffällig in Erscheinung tritt, liegt daran, 
daß eine eigentliche Umrechnung an keiner Stelle wirk- 
lich durchgeführt wird, daß es vielmehr bei der Emp- 
fehlung bleibt. Man vermißt auch Bemerkungen 
darüber, wo die Entsprechungen ihre Grenze finden, 
daß es mechanische Gebilde gibt (und welche), denen 
keine elektrischen Schaltungen entsprechen, und elek- 
trische Schaltungen, die sich nicht ins Mechanische ab- 
bilden lassen. 

Ganz zum Schluß sei noch bemerkt, daß die Er- 
örterungen über den Phasensprung bei einfachen Schwe- 
bungen den Verfasser zu einem Fehlschluß verleitet 
haben (siehe die „Kleine Mitteilung‘ in einem der 
nächsten Hefte). 

Den hier angemerkten Mängeln steht aber eine Fülle 
von ausgezeichneten Leistungen gegenüber, so daß.das 
Buch ganz zweifellos dennoch zu den besten auf diesem 
Gebiete gezählt werden muß. 


Karlsruhe, K. Klotter. 


Dr. phil. Georg Hamel (o. Prof. an der Technischen 
Hochschule Berlin, Integralgleichungen, 
Einführung in Lehre und Gebrauch. 
(2. beriehtigte Auflage.) VIII u. 166 S. mit 19 Abb. 
Berlin-Göttingen-Heidelberg 1949. Springer-Verlag. 
Preis brosch. 15,60 DM. 

Hier erscheint als nur wenig veränderter Neudruck 
ein Buch, das sichschon mit seiner ersten Auflage einen 
großen Leserkreis gesichert hat. (Eine ausführliche 
Besprechung steht in der ZAMM, Bd. 18 (1938) S. 204 
bis 205.) Es ist besonders dazu bestimmt, den Prak- 
tiker — Ingenieur und angewandten Mathematiker — 
in das wichtige Gebiet der Integralgleichung einzu- 
führen, das Interesse für diese elegante Theorie zu 
wecken und das Studium der Originalliteratur vor- 
zubereiten. Ohne große Vorkenntnisse wird im engen 
Zusammenhang mit praktischen Beispielen die lineare 
Integralgleichungstheorie von E. Schmidt und die 
von Fredholm recht weit durchgeführt. Daneben 
werden, dem oben genannten Ziel des Buches angepaßt, 
die Grundgedanken anderer Arbeiten, z. B. von Hil- 
bert und Enskog, wenigstens skizziert. 

Bestimmt wird auch die neue Auflage dieses leicht 
verständlichen, an keiner Stelle langatmigen oder er- 
müdenden Buches zahlreiche Studenten und Praktiker 
in das Gebiet der’ Integralgleichung einführen. 


N.J.Lehmann.. 


Jos. E. Hofmann, Die Entwicklungs- 
geschichte der Leibnizschen Mathe- 
matik während des Aufenthaltes in 
Paris (1672—1676). VIII u. 2528. mit 27 Abb. 
München 1949. Leibniz-Verlag. Preis geb. 26,— DM. 

Die hier vorliegende Darstellung der Entwicklung 
der Leibnizschen Mathematik während der rund 
dreieinhalb Jahre seines Aufenthaltes in Paris stützt 
sich auf reiches neues Material. Als Anfänger kam 
Leibniz dorthin. Die Pariser Jahre sind die Zeit, in 
der er seine entscheidenden Gedanken konzipierte, den 
Caleulus erfand, seine erste Rechenmaschine konstru- 
ierte, in deren letztem Jahre sich aber auch schon die 
Anfänge des Prioritätsstreites mit Newton abzeichnen. 
Es war eine äußerst schöpferische Periode, und der 
Verf. tut recht, sich in seiner Darstellung auf diese Zeit: 


Buch wird ein spannendes Bild der damaligen Ent 
lung der Mathematik und der daran beteiligten Pe 


a A 


Zusammenhänge, des Werdens und Entstehens der 
Kernprobleme, ihrer Gefüge, ihrer Wirksamkeit und 
Zielsetzung‘“ gegeben. Klar zeigt sich, daß hinsichtlich 
der Infinitesimalrechnung kein Plagiat von Leibniz 
vorliegt. Nachdem er, zunächst von Huygens be- 
raten, durch das Studium der ganzen ihm zur Verfü- 
gung stehenden Literatur die gefühlsmäßig richtige 
Kinstellung zu den Infinitesimalproblemen gewonnen 
hatte, hat er seine Methoden durchaus selbständig ent- 
wickelt. Dabei war das Entscheidende die Einführung 
einer zweckmäßigen Symbolik. 

Die Darstellung des Buches ist außerordentlich sorg- 
fältig. Jede Feststellung ist durch Hinweis auf Brief- 
stellen oder Stellen in den Leibnizschen Aufzeich- 
nungen usw. belegt, wenn diese 'Stellen selbst hier auch 
nicht abgedruckt sind. Ein chronologisches Verzeich- 
nis der rund 700 benutzten Briefe und Aufzeichnungen, 

. ein umfangreiches Namen- und Schriftenverzeichnis, 
sowie ein Sachverzeichnis findet man am Schluß dieses 
ausgezeichneten Buches. 


Dresden. \Wrllers, 

Dr.-Ing. Johannes Johannson, Das Cross-Ver- 
fahren. Die Berechnung biegefester 
Tragwerke nach der Methode des Mo- 
mentenausgleichs. 1238. mit 137 Abb. und 
18 Zahlenbeispielen. Berlin/Göttingen/Heidelberg1948. 

' Springer-Verlag. Preis 14,40 DM. 

Der Spannungszustand der statisch unbestimmten 
Stabwerke kann entweder mit statisch unbestimmten 
Schnittkräften unmittelbar entwickelt oder aus dem 
Verschiebungszustand des Tragwerks mit Hilfe von 
Knoten- und Stabdrehwinkeln abgeleitet werden. 
In beiden Fällen entstehen zunächst Ansätze von 
linearen Gleichungen, deren Auflösung oft durch die 
Fehlerempfindlichkeit Schwierigkeiten bereitet. Man 


'Buchbesprechungen i 
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Erleichterung bedeutet. Berücksichtigt sind insbeson- 
dere folgende Einflüsse und Varianten: Spannung, 
Dehnung, Eigengewicht, äußere Belastung, Wärme- 
dehnung, wechselnde Temperatur bei konstanter Be- 
lastung, wechselnde Belastung bei gleicher Tempera- 
tur, Durchhang. Voraussetzung für die vereinfachte 
Rechnungsweise, wie sie hier zugrunde liegt, ist die 
Gültigkeit des linearen (H ookeschen) Elastizitäts- 
gesetzes und die für nur geringen Durchgang geltende 
Approximation der Kettenlinie als Parabel. 
Dresden. H. Neu ber. 


Prof. Dr. Werner Schaub, Weltraumflug (astro- 
nomische und physikalische Grundlagen). 938. m. 
20 Abb. Bonn 1949. Ferd. Dümmlers Verlag. Preis 
brosch. 6,35 DM. 

Der zu allen Zeiten diskutierte und immer wieder 
reizvolle Gedanke der Möglichkeit des Weltraumfluges 
wird im vorliegenden Büchlein im wesentlichen vom 
reinen Standpunkt der Dynamik materieller Systeme 
behandelt. Gegenstand der Erörterungen des Ver- 
fassers sind daher vor allem die bekannten Folgerungen 


‚aus den Gleichgewichtsbedingungen einer Punktmasse, 


die sich im Gravitationsfeld der Sonne und ihrer Pla- 
neten entsprechend dem Trägheitsgesetz der klassi- 
schen Mechanik bewegt. Weitere Eınflüsse, wie z. B. 
der Luftwiderstand in Erdnähe sind außer Betracht 
gelassen. Eskam dem Verfasser offenbar nur darauf an, 
die für ein Raumschiff folgenden punktmechanischen 
Gesetzmäßigkeiten möglichst anschaulich und ohne 
mathematische Ableitungen einem größerem Leserkreis 
nahezubringen. Dieses Ziel Kann im wesentlichen als 
erreicht angesehen werden. Der mathematisch etwas 
mehr geschulte Leser findet überdies im Anhang einige 
mathematische Ergänzungen. Allerdings kann natür- 
lich die Frage nach der praktischen Verwirklichung des ° 


ice hm ern, Lane drch nen Welle rt eingehender Konad De 
BLCH N near H Oro A ndt Grinter@Z ueksichtigung der weiteren mannigfachen Neben- 


vorgeschlagen und später in der deutschen Literatur 
wiederholt behandelt worden. H. Cross verzichtet auf 
den üblichen Ansatz zur Berechnung der Knoten- und 
Stabdrehwinkel, gewinnt die Lösung vielmehr durch 
unmittelbare Iteration mit Hilfe der Abklingung der 
Stabendmomente. Die Lösung ist daher auch unter 
dem Namen Momentenausgleichsverfahren bekannt 
geworden. 
gehend beschrieben, die inneren Zusammenhänge un- 
tersucht und die Anwendung durch gut gewählte 
Zahlenbeispiele ausführlich erläutert. Dies geschieht 
in vier Abschnitten, welche die Stabwerke mit un- 
verschieblicher und verschieblicher Knotenpunkt- 
figur, die Einflußlinien der Schnittkräfte und die An- 
wendung der Lösung bei veränderlichem Trägheits- 
moment der Rahmenstäbe behandeln. Damit ist auch 
der späteren Entwicklung des Verfahrens ausführlich 
_  Reehnung getragen. Der Verfasser hat sich mit dieser 
Arbeit den Dank der zahlreichen Bauingenieure ver- 
- dient, die sich mit der Berechnung biegefester Stab- 
werke ex officio zu beschäftigen haben. 


Dresden. K. Beyer. 
Prof. Wilhelm Staufer und Oberreichsbahnrat Wer- 


ner Splett,,‚Durchhang und Zugspannungen 
von Starkstrom-Freileitungen. IVu. 


1949. Verlag Richard Pflaum. Preis 14,— DM. 
Bei der Vorausberechnung von Hochspannungs- 
_  leitungen und bei der Überprüfung vorhandener Netze 
_ ergibtsich für den Sachbearbeiter die Schwierigkeit, die 
 Wechselbeziehungen zwischen den auftretenden. Vari- 
_  anten in möglichst klarer und nur erträglichen Rechen- 
aufwand erfordernder Form darzustellen. Die Ver- 
fasser haben sich daher aus praktischen Bedürfnissen 
ıs veranlaßt gesehen, für das Problem eihe Zu- 
lung von Kurventafeln und Tabellen zu 
h relche für die Praxis eine wesentliche 


J. Johannson hat das Verfahren ein- 


10 8., 32 Zahlentafeln u. 9 Kurventafeln. München 


' Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


bedingungen diskutiert werden, so daß das eigentliche 
Problem nach wie vor ungeklärt bleibt. 


Dresden: . H. Neuber. 


Dr. Fritz Chemnitius, Kleines Handbuch 
der Infinitesimalrechnung. VIII+184S. 
Pößneck 1949. Rudolf A. Lang-Verlag. Preis geb. 
14,— DM. 

Das Buch enthält auf seinen rechten Seiten 350 in 
allen Einzelheiten durchgerechnete Integrale, die größ- 


- tenteils auf elementare Funktionen führen. Auf den 


linken Seiten werden die erhaltenen Funktionen wieder 
differentiiert. So ist eine Sammlung häufig auftreten- 
der Integrale entstanden (Ref. konnte sie bereits mehr- 
fachim Lehrbetrieb verwenden), bei. der der Suchende 
nicht nur das Ergebnis, sondern auch die anzuwen- 
dende Methode findet. Allerdings ist der Zweck der 
Probe durch Differentiation des Ergebnisses nicht ganz ' 
einzusehen, da das Buch für ‚‚Leser mit einer gewissen 
Sachkenntnis‘ bestimmt ist, dem diese Rechnung keine 
Schwierigkeiten bereiten dürfte. Durch Fortfall dieses 
Teiles wäre Raum für weitere Integrale frei geworden. 
Atıßerdem hätte sich eine gewisse Straffung erreichen 
lassen durch die Auslassung einer großen Zahl von Inte- 
gralen, die Sonderfälle von anderen durchgerechneten 
Beispielen sind (z. B. rationale Funktionen von e“ als 
Sonderfälle der gleichen Funktionen von @). Bedenk- 
lich erscheint es, wenn bei Integration. durch Potenz- 
reihenentwicklung unterlassen wird, den Konvergenz- 
radius anzugeben. Aus dem Rahmen des Buches fallen 
16 willkürliche Beispiele für die Integration totaler 
Differentiale und die mehrmalige Integration einer 
Funktion über die gleiche Veränderliche, Im ganzen 
entsteht so der Eindruck, als wäre dieses an sich brauch- 
bare Buch mehr das Ergebnis einer Liebhaberei als 


‚systematischer Beschäftigung mit Infinitesimalrech- 


nung. 


Dresden. , G. Opitz. 
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Eingegangene Bücher — Zuschriften an den Herausgeber 


Z.angew. Math. Mech. 
Bd. 30 Nr.4 April 1950 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Optik-Photo und Kino-Feinmecha- 
nik - Uhren - Medizinmechanik und 
Elektromedizin 1949/50 (Warenverzeich- 
nis und Hersteller, Bundesrepublik Deutschland und 
West-Berlin). Herausgegeben in Zusammenarbeit mit 
dem Verband der Deutschen Feinmechanischen und 
Optischen Industrie. Wiesbaden 1949. Verlag für 
Wirtschaftsschrifttum Otto K. Kraußkopf. 


Dr. P. Luckey (Oberstud.-Rat i.R.), Nomogra- 
phie (Mathematisch-Physikalische Bibliothek Reihe I 
Nr. 59/60). 6. Aufl. 1078. m. 57Abb. Leipzig 
1949, Verlagsgesellschaft B. G. Teubner. Preis kart. 
4,20 DM. 


H. Boltz, Formeln und Tafeln zur nu- 
merischen (nicht logarithmischen) 
Bereehnung Gauß-Krügerscher Ko- 
ordinaten aus den geographischen 
Koordinaten, (Veröffentlichung des Geodäti- 
schen Institutes Potsdam. Neue Folge Nr. 11.) XVI 
+ 708. Potsdam 1943. 


F.Mühlig, Der 24m-Interferenzkompa- 
rator des Geodätischen Institutesin 
Potsdam. (Deutsche Akademie der Wissenschaften 
zu Berlin, Veröffentlichungen des Geodätischen Insti- 
tutesin Potsdam Nr. 2.) IV + 50 S.m.30 Abb. Berlin 
1949. Akademie-Verlag. 


W.Jenne, ZurAuflösunglinearerGlei. 
ehungssysteme. (Dautsche Akademie der Wis- 
senschaften zu Berlin. Geodätisches Institut Potsdam.) 
238. mit 25 Abb. Abdruck aus den Astronomischen 
Nachrichten. Bd. 278, Heft 1/2. Berlin 1949. 


Dr.-Ing. habil. Ernst Eckert, Einführung in 
den Wärme- und Stoffaustausch. VII 
+ 2038. m. 125 Abb. Berlin-Göttingen-Heidelberg 
1949. Springer-Verlag. Preis brosch. 21,— DM, geb. 
24,— DM. 


Dr.-Ing. habil. Walther Kauffmann (o. Prof. a. d. 
Techn. Hochschule in München), Einführungin 
die Teehnische Mechanik. Erster Band: 
Statik starrer Körper. VI + 166 S. m. 194 Abb. Ber- 
lin-Göttingen-Heidelberg 1949. Springer-Verlag. Preis 
brosch. 15,— DM. 


Prof. Dr.-Ing. A. Leonhard, Die selbsttätige 
Regelung (Theoretische Grundlagen mit prakti- 
schen Beispielen). IX + 284 S. mit 254 Abb. Berlin- 
Göttingen-Heidelberg 1949. Springer-Verlag. Preis 
brosch. 24,— DM, geb. 27,— DM. N 

Dr.-Ing. Winfried Oppelt, Stetige Regelvor- 
gänge (Bücher der Technik. Herausgeber: Dr.-Ing. 
Alfred Kuhlenkamp.) 1448. m. 42 Abb. Hannover- 
Wolfenbüttel 1949 Wissensch. Verlagsanstalt K. G. 
Hannover und Wolfenbütteler Verlagsanstalt GmbH. 
Preis brosch. 7,— DM, geb. 7,80 DM. 


Dr.-Ing. 6. Grüss (o. Prof. a. d. Bergakademie Frei- 
berg), Differential- und Integralrech- 
nung (Mathematik und ihre Anwendung in Physik 
und Technik. Herausgegeben von E, Kamke und 
A. Kratzer. Reihe A. Bd. 21). XVI + 642 S. mit 355 
Abb. Leipzig 1949. Akademische Verlagsgesellschaft. 
Preis geb. 39,— DM. 


Dr.-Ing. Th. Pöschl (o. Prof. a. d. Techn. Hochschule 
Karlsruhe), Lehrbuch der Technischen 
MechanikI: StatikundDynamik(3.um- 
gearbeitete Auflage.) VIIL + 343 S. m. 257 Abb. Ber- 
lin-Heidelberg-Göttingen 1949. Springer-Verlag. Preis 
geheftet 22,50 DM, geb. 25,— DM. 


Dr. K.L. Schecher und Studienrat Dr. Karl Seebach 
(Privatdozent a. d. Techn. Hochschule München), 
Einführung in die Mathematik, ein 
Lehrbuch für Selbstunterriehtung und Wiederholung. 
Bd.I. 288 S. m. 408 Abb. Frankfurt/M. 1950. Verlag 
Christian Schmidt. 


ZUSCHRIFTEN AN DEN HERAUSGEBER. 


Zu R. Zurmühl: Zur numerischen Auflösung 
linearer Gleichungssysteme nach dem Matri- 
zenverfahren von Banachiewicz. Z. angew. 
Math. Mech. 29 (1949) 76 — 84. 


In mehreren Arbeiten hat Herr Banachiewicz ein 
Verfahren zur Lösung linearer Gleichungen aufgestellt, 
das er als wesentlich besser ansieht als das Gaußsche 
Eliminationsverfahren; es beruht darauf, die Glei- 
chungsmatrixA in das Produkt zweier entgegengesetzter 
Dreiecksmatrizen zu zerlegen. Die Methode ist bisher 
kritisch nie betrachtet worden, so daß ich annahm, ihr 
Kern sei allgemeiner erkannt. So unterließ ich in einer 
früheren Arbeit über lineare Gleichungen!) eine Kritik. 

Zum ersten Male wurde ich durch eine Arbeit von 
Cassina?) darauf aufmerksam, daß man den Kern der 
Methode nicht erkannt hatte. Denn trotzdem nach 
Cassina die Zahl der Multiplikationen und Additionen 
bei Banachiewiez und Gauß dieselbe ist und er in 
der Einleitung auch behauptet, daß der ‚praktische 
Vorteil bei Banachiewicz darin besteht, nur die notwen- 
digen Zahlen der Elimination beizubehalten, ohne die 
Elimination wirklich auszuführen‘, bleibt er eine nähere 


%) Bericht über die verschiedenen Methoden zur Lösung eines 
Systemslinearer Gleichungen mitreellen Koeffizienten, Akademie 
VE aBnHaDpen Amsterdam, Proceedings, Bd.50 u, 51, 1947 u, 
1948, 

?) Numero delle operazioni elementari necessarie per la riso- 
luzione dei ssistemi di equazioni lineari, Boll, Un, Mat, Ital, (3) 3, 
1948), 3. 142—147, . 


Ausführung dieser Behauptung schuldig, beschreibt 
vielmehr die Methode in herkömmlicher Weise durch 
die Multiplikation zweier Dreiecksmatrizen. 

Auch in der Arbeit von Herrn Zurmühl wird der 
Zusammenhang nicht gesehen. Im Gegenteil wird aus- 
drücklich eine wesentliche Verschiedenheit beider Ver- 
fahren betont, ferner die große Wichtigkeit des Ver- 
fahrens von Banachiewicz hervorgehoben, das nur der 
Krakovianen-Schreibweise seine Unbekanntheit ver- 
danke, undschließlich demVerfahren eine entscheidende 
Herabsetzung des Arbeitsaufwandes zugeschrieben. 


Diese Irrtümer erklären sich dadurch, daß die Zer- 
legung der Matrix A in das erwähnte Produkt, also 
DJ=NN, mit Elimination nichts zu tun zu haben 
scheint. In Wirklichkeit überlegt man sich aber leicht, 
daß der sich aus der Produktdarstellung von A er- 
gebende Algorithmus genau der gleiche ist wie der des 
sog. abgekürzten Gaußschen Verfahrens. Wenn 
man nämlich von jedem abgeleiteten Gleichungssystem 
stets eine feste, etwa die erste, Gleichung herausgreifen 
will, so braucht man nicht von jedem System alle 
Gleichungen zu bestimmen und aufzuschreiben, viel- 
mehr genügt dazu jedesmal die erste Gleichung und 
allerdings noch die erste Spalte. Es ist wohl unnötig, 
diesen einfachen Algorithmus nochmals zu beschreiben. 
Er hat — obwohl mit dem von Banachiewicez völlig 


identisch — den Vorteil, daß er infolge seiner Herkunft 2 


nie vergessen werden kann. In der Arbeit von Herrn 
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Zurmühl ist also in dem Schema auf $. 82 die erste 
Zeile(abgesehenvomLeitkoeffizienten)gleich der ersten 
Gleichung des ersten Systems und ist die erste Spalte 
gleich der ersten Spalte des ersten Systems. Ferner ist 
die zweite Zeile rechts von der Diagonalen gleich der 
ersten Gleichung des ersten abgeleiteten Systems und 
ist die zweite Spalte unterhalb der Diagonalen gleich 
der ersten Spalte des ersten abgeleiteten Systems usw. 

Diese Ableitung hat aber den weiteren Vorteil, daß 
man Störungen beim Funktionieren des Algorithmus 
ohne weiteres ausschalten kann, während bei Bana- 
chiewicz dazu extra Betrachtungen nötig sind. Wenn 
nämlich der Leitkoeffizient eines abgeleiteten Systems 
zu klein wird (oder verschwindet), so kann man ent- 
weder die Variablen umstellen oder eine andere Glei- 


- chung anstatt der ersten nehmen. 


> 


Die Zerlegung A=CBführtalsonurscheinbar 
zu einer andern Auflösungsart als die Elimi- 
nationsmethode. 

Aber sogar die beiden Matrizen C und B werden bei 
dem E iminationsverfahren gleichzeitig mitgeliefert. 
Bei Matrix B ist dies klar, denn sie ist die Matrix des 
dreieckigen Endsystems, und die Spalten von Matrix C 
erhält man nach Obigem als die ersten Spalten des 
abgeleiteten Systems. Ohne Bezugnahme jedoch auf 
den aus A=CB folgenden Algorithmus ergibt sich der 
Beweis dafür folgendermaßen. 

Mit Hilfe von Matrizen drückt sich der Eliminations 
prozeß bekanntlich!) so aus, daß 


Rn-ıEn-a''' RR,A=L ) 2eilehlte 'e (1), 


wo L und die R;linke bzw. rechte Dreiecksmatrizen sind 
(mit Nullen links bzw. rechts der Diagonalen) und z. B. 


120.020 0 

—.0944,.10.0 0 
RK=|-a, 01 0813 

N 

1 DR, Sl) 

0 N) 0 
R=|0 —ay 1 0% 


z Ba BES 
Der Fehler bestand nun bisher darin, zur Zerlegung von 


A= RTL zuerst das Produkt Ra-ı *‘* Rekı=R 
zu bestimmen und dann die Reziproke R-!, während 


man in Wirklichkeit erst die Reziproken R;® zu bilden 


hat und dann deren Produkt. Zerlegt man nämlich 
jedes R; in seinen Diagonalteil Z und den außerhalb 
davon liegenden Teil U;: 


ne BU Ne (2% 

so findet man sofort, daß 
U} =0= Nullmatrix 
U;U;=0=Nullmatrix für k>i.. ..(4). 


Daher ist 
" BE=E—U;=(E+UVi)(E—Uı), 


also 


Ri =(ELU)"=E—UVi....(d), 


‚daher nach (l):. 


a 
ASRRZ Rn 
= (B—=U,)(E—U,) -- (E—-Unn)L 


| und nach (4): 


i N a a re N TB 


oder, da die einzige Spalte von — U; gleich der ersten 


Spalte si des i-ten Gleichungssystems ist: 


A=(E+as+s+ ++ 81)ZL. 


ERRE 


a8. 8. 0.2), 8. 937/38. 
Ca 7 TR n a 


N 


Die Klammer ist dabei nichts anderes als die Matrix C 
in der Bezeichnung von Herrn Zurmühl. 


Ergebnis. 


1. Das Verfahren von Banachiewicz ist kein selb- . 


ständiges Verfahren, sondern ist identisch mit dem 
abgekürzten Verfahren von Gauß. 
2. Einer Zerlegung A=CB bedarf es dazu nicht 


ausdrücklich. Doch werden beide Matrizen B, C auch 


im Gaußschen Algorithmus automatisch mitgeliefert. 
Sie bestehen aus der jeweils ersten Zeile bzw. ersten 
Spalte der einzelnen abgeleiteten Systeme. 

3. Die Zerlegung A=CB zum Ausgangspunkt der 
Betrachtungen zu nehmen, ist sonderbarerweise sogar 
hinderlich für ein Erkennen des Mechanismus der 
Lösung. 

4. Auch die Methode von Choleski ist demgemäß 
eine Variante der Gaußschen Eliminationsmethode. Sie 
ergibt sichdurch andereWahl der Diagonalkoeffizienten 
von C und B, nämlich durch Gleichmachen entspre- 
chender Elemente. Auch diese Methode wird daher 
besser und klarer als Eliminationsmethode durch- 
geführt. 


den Haag. E. Bodewig. 


Erwiderung 

Das von mir in Bd. 29, S. 76—84 dieser Zeitschrift 
beschriebene Matrizenverfahren von Banachie- 
wicz zur Auflösung linearer Gleichungssysteme ist 
den zugrunde liegenden Formeln nach in der Tat 
identisch mit dem sogenannten abgekürzten Gauß- 
schen Verfahren. Die gleichen Überlegungen wie 
Herr Bodewig führt H. Unger, Darmstadt, 
in einer noch unveröffentlichten Arbeit durch. Auch 
sonst: findet: sich ähnliches schon mehr oder weniger 
klar ausgesprochen!). Dagegen\enthält der in obiger 
Zuschrift Anm.1 zitierte Bericht von Bodewig bei 
der Besprechung des gleichartigen Verfahrens von 
Cholesky (1947, S.939) noch nichts darüber, und 
erstin einem Nachtrag (1948,S. 214) wird das Ergebnis 
vonAndersen angeführt. Daß der Zusammenhang 
zwischen dem Gaußschen Algorithmus und der 
neuen Rechenanordnung nunmehr in so vollständiger 
Weise aufgedeckt ist, ist sehr zu begrüßen. Um so 
vertrauenswürdiger erseheint damit auch das neue 
Vorgehen. — Im Falle kleiner oder verschwindender 
Leitkoeffizienten hilft man sich beim Matrizenver- 
fahren nicht anders als bei Gauß, nämlich durch 
Umstellen der Variablen oder der Gleichungen, vgl. 
Nr. 8 meiner Arbeit. 

Obgleich nun das abgekürzte Gaußsche Ver- 
fahren seit langem bekannt ist, so ist doch die in ihm 
gelegene Möglichkeit einer weiteren Rechnungsverein- 


. fachung durch volle Ausnutzung der Rechenmaschine 


lange ganz unbeachtet geblieben. Keines der ein- 
schlägigen Lehrbücher macht davon Gebrauch. Erst 
der auch hier als heuristisches Prinzip sich bewährende 
Matrizenkalkül läßt diese Möglichkeit ins Bewußtsein 
treten, zuerst bi Cholesky, später dann bei 
Banachiewicz. Daß sich bei beiden nachträg- 
lich die engste Verwandtschaft mit dm Gauß- 
Algorithmus herausstellt, kann der großen praktischen 
Bedeutung der neuen Vorgehensweise keinen Ab- 
bruch tun.. Denn das Problem der Gleichungsauf- 
lösung ist, wenn man von den Iterationsverfahren 
absieht, gar kein theoretisches, sondern seinem inner- 
sten Wesen nach ein rein praktisches. Jede Variante, 


welche den effektiven Arbeitsaufwand entscheidend 
- herabsetzt, muß als ein Fortschritt gewertet werden, 


auch dann, wenn die ihr zugrunde liegenden Formeln 
nicht neu sind. 

Es handelt sich hier auch nicht um Trivialitäten, 
Gewiß ist: der Schritt vom abgekürzten Gauß- 


ı) H, Tensen, An attempt at a systematie classification of 
some methods for the solution of normal equations, Kopenhagen 
1944, — E. Andersen, Solution of great systems of rormal 
equations together with an investiration of Andraes dot-figure, 
An arithmetical-technical investigation. Kopenhagen 1947. 
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schen Verfahren bis zur neuen Rechenanordnung 


nicht groß; doch ist er, wie die historische Entwick- 
lung gezeigt hat, auch nicht selbstverständlich. Die 
große ihm zukommende Bedeutung findet sich im 
Schrifttum nur selten deutlich ausgesprochen, und noch 
der oben zitierte Bericht von Bode wig kommtdurch 
das bloße Abzählen der Operationen zu einer unrich- 
tigen Einschätzung des Verfahrens von Cholesky. 
Auf die neuartige Rechenanordnung und die sich aus 
ihr ergebende außerordentliche Zeitersparnis mit 
Nachdruck hinzuweisen, war das Ziel meiner Arbeit. 
Die dort behandelten Gesichtspunkte behalten auch 
nach Aufhellung der theoretischen Zusammenhänge 
ihre unverminderte Gültigkeit. 


Darmstadt. R. Zurmühl. 


Zu G. Opitz: Praktische Verfahren zur Lösung: 
von Gleichungen vierten Grades. Z. angew. 
Math. Mech. 25/27 (1947), S. 171/173 und zu 
H. Blenk: Nomogramme für die Gleichung 
4. Grades mit reellen oder komplexen Wurzeln. 
Z. angew. Math. Mech. 29 (1949), S. 58/61. 


Zu den von G. Opitz angegebenen Verfahren sei 
eine Variante mitgeteilt. Diese geht ebenfalls davon 
aus, die Gleichung 


J[e))=z"+,. +, ? +2 +9,=0 
in zwei quadratische Faktoren zu zerlegen: 
[)=(—s12+ 9): ( — 324 9) =0. 
Die Größen s}, ?1, 82, 9, sind stets reell. Für sie 
— wie auch Opitz angibt: 
PPR=-% 178=—4, 
PtrPm)ts32 =, PRA4tmMp=—G. 
Die Variante, die mir vor Jahren unabhängig von- 
einander von den Göttinger Herren Dr. Wielandt 
und Dr. Koch mitgeteilt wurde, besteht nun darin, 
daß man verschiedene Werte p, aus dem Bereich 
-Nal<m<+ Val 


annimmt und dafür jeweils ermittelt: 


a Ag — 4A 

LG 3 121 

PP; en = is 
PL pı Pa 


PH +Pm)+ 3% =A. 


Derjenige Wert p, ist dann der richtige, für den 
A=.a, ist. Durch systematisches Probieren (Ein- 
gabeln) kann man diesen Wert mit jeder gewünschten 
Genauigkeit bestimmen. 


‚ Die zu diesem Werte p, gehörigen Werte p,, 1, 8, 
liefern dann — durch Auflösen der beiden quadra- 
tischen Gleichungen: 


2 — 2 +9p=0 und 2—-,2+9%,=0 


die vier gesuchten Wurzeln von f(z) = 0. 


Dabei möchte ich zur Auflösung dieser quadratischen 


Gleichungen auf eine Formel hinweisen, die nicht 
allgemein bekannt zu sein scheint. Diese Formel, die 
mir früher Herr Warning, Hamburg, mitgeteilt 
hat, umgeht die Fehlerempfindlichkeit, die die nor- 
male Auflösungsformel mitunter aufweist. Danach 
hat: die Gleichung 


ha, =0 
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St. Monino bei Moskau. KarlBorkmann. 


Kritische Betrachtungen zur sogenannten 
Großzahlforschung. 


Uber die Arbeit von P. Riebesellin Band 28, 
Heft 7/8 dieser Zeitschrift (1948) hat zwischen Herrn 
Riebesell und mir eine längere Diskussion stattgefun- . 
den, deren Ergebnis ich in einem Schreiben an Herrn 
Riebesell folgendermaßen zusammengefaßt habe: 

1. Ihr Artikel hatte nicht den Sinn, eine technische 
Beurteilung statistischer Daten abzulehnen, sondern 
vor allem darauf hinzuweisen — was bislang meines 
Wissens in der deutschen Literatur zu wenig beachtet 
wurde —, daß ein statistischer Schluß ohne weiter- 
gehende Untersuchung gefährlich ist und ein derartiges 
Urteil besser nicht als statistische Schlußfolgerung, 
sondern ehrlich als technisches Urteil bezeichnet werden 
sollte, auch wenn bei seiner Bildung statistische Me- 
thoden mitgenutzt wurden. s 

2. Bezüglich der Freiheitsgrade ... auditatur et 
altera pars, ich fürchte, jede Jury wird einen dissenting 
judge aufweisen. 

3. Bezüglich der Bewertung ‚genügend kennzeich- 
nend‘‘ werden die Meinungen immer auseinandergehen. 
Geben zwei Bearbeiter auf Grund der subjektiven Be- 
wertung der Folgen eines Fehlurteils das Urteil auf 
Grund verschiedener &- oder o-Grenzen, so kann einer 
sagen „genügend kennzeichnend‘‘ der andere „nicht- 
kennzeichnend‘‘, ohne daß einem ein statistischer Vor- 
wurf gemacht werden könnte. Schematische Fest- 
legung einer Grenze — sei es die von mir für Erstent- 
scheidungen bevorzugte & = 0,05 bzw. 0,95 Grenze, 
sei es die klassische 3 -Grenze — dürfte allerdings in 
keinem Fall zweckmäßig sein. 

4. In gleicher Weise gelten technische Erwägungen 
für den Ansatz kennzeichnender Tests. Je nach Be- 
wertung des Einzelfalles wird bewußt oder unbewußt 
auf eine vielleicht mögliche Trennung in Teilkollektive 
verzichtet werden, wenn ein globaler Test wie der y?- oder 
Mutungsgrenzansatz das Vorliegen eines einheitlichen 
Kollektives nicht kennzeichnend unwahrscheinlich 
macht und schärferes Urteil für den Untersuchungs- 
zweck nicht erforderlich ist. 


Geben schärfere Tests für zwei verschiedene Null. 
hypothesen unterschiedliche Wahrscheinlichkeiten, die 
aber alle nicht kennzeichnend sind wie ım Zahlenbeispiel 
die 30-Grenze, so ist: eine oft rein technisch bedingte 
Wahl erforderlich, wenn die Aufgabenstellung die 
größere Aussageschärfe fordert: Unterschiede im & 
innerhalb des als nicht kennzeichnend erachteten Be- 
reichs berechtigten weder zur Annahme noch zur Ab- 
lehnung einer Hypothese. Wird die größere Aussagen- 
schärfe nicht: verlangt oder in Zukunft nicht nach- 
geprüft, so ist auch nach meiner Auffassung Verzieht/n. | 
auf Trennung und Beschränkung auf die geringt 
mögliche Parameterzahl am Platze: die techni she 
Aussagenschärfe kann zwar leiden, die Sicherheit sta 
stistischer Voraussagen gewinnt aber oft. 
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